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Abstract 

According to Belinskii, Khalatnikov and Lifshitz (BKL), a generic spacelike singularity 
is characterized by asymptotic locality: Asymptotically, toward the singularity, each spatial 
point evolves independently from its neighbors, in an oscillatory manner that is represented by 
a sequence of Bianchi type I and II vacuum models. Recent investigations support a modified 
conjecture: The formation of spatial structures ('spikes') breaks asymptotic locality. The 
complete description of a generic spacelike singularity involves spike oscillations, which are 
described by sequences of Bianchi type I and certain inhomogeneous vacuum models. In this 
paper we describe how BKL and spike oscillations arise from concatenations of exact solutions 
in a Hubble-normalized state space setting, suggesting the existence of hidden symmetries and 
showing that the results of BKL are part of a greater picture. 



1 Introduction 

This paper is part of a research program on generic spacetime singularities. It is concerned with 
spike oscillations, a missing piece in the description of the asymptotic dynamics of spacetimes 
toward a generic spacelike singularity. 

In research beginning in the 1960's, Belinskii, Khalatnikov and Lifshitz (BKL) claimed to have con- 
structed an approximate general solution to Einstein's field equations in the vicinity of a spacelikc 
('cosmological') singularity [1, 2, 3]. The central assumption of BKL is that the asymptotic dy- 
namics is 'local' in the sense that each spatial point evolves toward the singularity individually and 
independently of its neighbors as a spatially homogeneous model. In addition, BKL gave heuristic 
arguments that suggested that perfect fluid models with soft equations of state, such as dust or 
radiation, are asymptotically 'vacuum dominated', i.e., asymptotically toward the singularity the 
spacetime geometry is not influenced by the matter content. Moreover, the BKL conjecture stated 
that for generic asymptotically vacuum dominated models the evolution of each spatial point in the 
vicinity of the singularity is approximated by a sequence of vacuum Bianchi type I solutions (i.e., 
Kasncr solutions) mediated by vacuum Bianchi type II solutions, where the latter determine a dis- 
crete 'Kasner map' with chaotic properties [4]-[21]. In general, the Kasner map generates infinite 
sequences of Kasner states, and BKL argued that these sequences describe oscillatory temporal 
behavior toward a generic spacelike singularity. 

In 1992 Belinskii pointed out that the independent evolution of different timelines (due to the 
locality assumption) leads to a lack of temporal synchronization of oscillations, and hence spatial 
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structure formation, reminiscent of turbulence, takes place, which has raised the issue if the locality 
assumption of BKL is consistent [22], see also [19] and references therein. 

In 1993 Berger and Moncrief numerically studied vacuum G2 models (i.e., models with two space- 
like commuting Killing vectors [7]) with T 3 topology, so-called T 3 Gowdy models [23, 24], and 
observed the development of large spatial derivatives near the singularity, which they called 'spiky 
features' [25] . Furthermore, based on work by Grubisic and Moncrief of the same year [26] , these 
features were tied to a condition at isolated spatial surfaces that could be imposed on a formal 
expansion of the Gowdy metric near the singularity. 1 

Toward the end of the 1990's, Berger, Moncrief and coworkers had found further numerical evidence 
that the BKL picture seemed to be correct generically in vacuum G2 spacetimes, and even in 
spacetimes with only one spacelike Killing vector, but that there were difficulties at exceptional 
isolated spatial surfaces [27]-[30], see also [10] and references therein. About the same time, Hern 
independently resolved individual spatially spiky features to high numerical accuracy, although for 
short time scales [31, 32]. 

A significant analytic step toward an understanding of spiky features in T 3 Gowdy models was 
accomplished in 2001 by Rendall and Weaver [33]. They applied a solution generating technique 
and Fuchsian methods in [34, 35], to produce asymptotic expansions for 'spikes', which they referred 
to as 'true' or 'false spikes', where false spikes where shown to be gauge artefacts. The work on 
spiky features in Gowdy spacetimes by Rendall and Weaver was followed up by Garfinkle and 
Weaver in 2003 who used two different complementary numerical techniques [36]. In particular 
they studied so-called 'high velocity spikes' and found that they eventually transform into 'low 
velocity spikes'. 2 

Up to this point, the framework for studying singularities in inhomogcncous spacetimes had been 
dominated by the metric approach, as exemplified by the work of BKL, and the Hamiltonian 
approach used by Moncrief and coworkers. However, at the same time, the studies of spatially 
homogeneous models were dominated by the 'Hubble-normalizcd' dynamical systems approach 
introduced by Wainwright and Hsu [38], which resulted in substantial progress [7, 39, 40]. 3 In 
2003 in an attempt to generalize this framework to a general inhomogeneous context, Uggla et 
al [46] (UEWE) recast Einstein's field equations into an infinite dimensional dynamical system by 
means of 'Hubble-normalized' scale-invariant variables. This resulted in a more specific and precise 
formulation of the BKL conjecture in terms of the asymptotically regularized field equations on 
an asymptotically bounded state space. The approach in UEWE was subsequently reformulated 
somewhat by means of a conformal transformation [47] which led to a geometric framework that 
was in turn specialized by means of a conformally Hubble-normalized Iwasawa frame by Heinzlc 
et al [48] to provide a link to work on cosmological billiards [49], and to establish the consistency 
of the BKL picture by means of dynamical systems methods. 

However, billiard and BKL consistency do not exclude other behavior at generic singularities. Nu- 
merical studies based on the dynamical systems formulation in UEWE gave support for the BKL 
picture generically [50, 51], and provided evidence that the structure formation predicted by Be- 
linskii [22] is not a threat to the consistency of the BKL scenario, since it develops on superhorizon 
scale toward the singularity (see also [52]). On the other hand, numerical and heuristic arguments 
also suggested that there exist recurring temporal spike oscillations along certain timelines. These 
timelines correspond to spatial points that form two-dimensional ' spike surfaces' 1 and their shrink- 
ing neighborhoods. The temporal spike oscillations are gauge invariant features that are distinct 
from BKL oscillations, as illustrated by the evolution of the Weyl scalars; in contrast to the BKL 
case, spatial derivatives contribute significantly to the Weyl scalars along the timelines of the spike 

x Due to the symmetries of G2 models, fields only depend on time and a single spatial coordinate; hence, two- 
dimensional symmetry surfaces, defined by a fixed value of this spatial coordinate, are often referred to as spatial 
points in the G2 literature. However, since our ultimate objective is generic singularities (in spacetimes without 
symmetries), and because in that context this nomenclature is inappropriate, we will refrain from referring to 
surfaces as points. 

2 Velocities and other diagnostic tools for describing and understanding asymptotic features associated with 
singularities are discussed in Appendix A. 4. For further and recent primarily analytic work on T 3 Gowdy models 
we refer to the review [37] and references therein. 

3 For recent developments about generic singularities in spatially homogeneous models, see [41]-[45]. 
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surfaces (and in their shrinking neighborhoods), and in this sense oscillatory spike evolution is 
'non-local'. 

Based on the solution generating technique used by Rendall and Weaver to produce approximate 
solutions [33] , in 2008 Lim was able to find a 1-parameter family of explicit (exact) inhomogeneous 
vacuum G2 1 spike solutions', which are expressible in terms of elementary functions [53]. As 
suggested by the work in [50, 51], these solutions form the building blocks for spike oscillations. 
This was further substantiated by Lim et al in 2009 [54] where the analytic spike solutions found 
by Lim were numerically matched with high accuracy to actual G2 solutions in the regime toward 
the singularity (at least for a small number of spike oscillations) . Furthermore, Lim et al managed 
to provide numerical evidence that showed that 'higher order spike solutions,' with several spike 
surfaces, could be described, locally and asymptotically, by the simpler 'first order spike solutions' 
with a single spike surface. Hence the first order spike solutions seem to be the building blocks for 
describing spike oscillations in general (at least for the spiky non-BKL behavior discovered so far 
in numerical simulations), and it is therefore these solutions that we will focus on in this paper. 4 

In the BKL scenario, Bianchi type II vacuum solutions ('transitions') connect Kasner solutions, 
and thus generate the oscillatory BKL behavior through heteroclinic chains (BKL chains). A 
specific goal of the present paper is to show how the explicit spike solutions by Lim connect the 
(generalized) Kasner metrics in the Hubble-normalized state space picture. Thereby we obtain a 
networks of 'spike chains' that represent the asymptotic dynamics of solutions and thus the spike 
oscillations. Throughout this paper we focus on general vacuum G2 models. There are reasons 
to believe that these models describe essential features of generic spacelike singularities, i.e., the 
general singularities that occur in models without symmetries [50, 54]; we will return to this issue 
in the concluding remarks. 

The paper is organized as follows. Section 2 provides a succinct overview of the G2 models in the 
conformal orthonormal (Iwasawa) frame approach. Section 3 is devoted to the essential building 
blocks of local BKL and non-local spike dynamics: Transitions. We discuss frame and curvature 
transitions and give a detailed account of the explicit high and low velocity spike solutions found by 
Lim tailored to our purposes and in the conformally Hubble-normalized picture. Section 4 treats 
a special class of G2 models, the orthogonally transitive G2 models which comprises, e.g., the T 3 
Gowdy models. We describe in detail how transitions can be concatenated to obtain (finite) chains 
of transitions that are relevant for the asymptotic dynamics of solutions. Finally, section 5 treats 
the general G2 case. In this case, concatenation leads to (in general) infinite BKL chains and spike 
chains that determine the oscillatory behavior of G2 solutions toward the singularity. Section 6 
contains the conclusions together with a brief description and assessment of previous numerical 
results in view of our present analytical and numerical understanding. We also discuss connections 
between BKL and non-BKL behavior and give hints to underlying hidden asymptotic symmetries. 
The appendices contain a discussion of on topological issues and a brief review of the results on T 3 
Gowdy spacetimes; we stress that the 'permanent' spiky features characteristic of the T 3 Gowdy 
case are replaced by transient recurrent features in the general G2 case: Spike oscillations. 

2 G2 models in the conformal orthonormal frame approach 

In this paper we study vacuum G2 models, which are characterized by the existence of two spacelike 
commuting Killing vector fields. In symmetry-adapted local coordinates these Killing vectors are 
d x and d v , which allows the metric to be written as 

ds 2 = -N 2 (dx ) 2 + e- 2bl (dx + mdy + n 2 dzf + (dy + n 3 dz) 2 + er 2b * dz 2 , (1) 

where the appearing functions depend on x° and the spatial coordinate z alone, cf. [56]. With 
n.i = — ni, n 2 = — n 2 + ^1^3, "3 = — n 3 we obtain the alternative form [30, 57, 58, 59] 

ds 2 = -N 2 (dx°f + er 2bl (dx - n x dy + (-n 2 + n^) dzf + e" 26 " (dy - n 3 dz) 2 + e~ 2b3 dz 2 . (f) 
4 For the influence of spike solutions on matter, see [55]. 
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We refer to Appendix A.l for details and for a topological discussion. The form (1) provides an 
'Iwasawa' representation of the metric and is naturally associated with the orthonormal spatial 
(co-)frame {ui , lo 2 , u> 3 } determined by 

e- bl 

e-» 2 | | ! n :i | j j . (2a) 


which is the dual of the orthonormal spatial 'Iwasawa frame' {ei, e 2 , &z\ (see, e.g., sections 2 and 13 
in [48]) given by 

/e fcl \ /l ni n 2 \ /cV 

V 
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e b I I 1 n 3 I \dy\. (2b) 




Combining this spatial frame with eo = N~ 1 dt and w° = Ndx° yields an Iwasawa based orthonor- 
mal frame for the spacetime metric, i.e., ds 2 — — u° <g> cj° + 5 a p io a <g> ujP , where we use Greek letters 
as spatial frame indices. 

The area density of the symmetry surfaces is R — er^^ . Whenever the gradient of R is timelike, 
it is possible to define time by choosing R to be a prescribed monotone function, e.g., R = x°, or 

-logR = b 1 + b 2 =t + t Q , (3) 

for some constant to (which is chosen to be to = unless stated otherwise) . As shown in [57, 60] , 
for the T 3 case, R oc exp(— t) G (0, oo) represents the maximal globally hyperbolic development 
of initial data at t = const, where t — > oo represents the initial singularity; the only exception is 
the Taub (flat Kasner) solution. In the area time gauge (3), in the nomenclature of [61], the line 
element can be written as 

ds 2 = e- 2b \-ae~ 2t dt 2 + dz 2 ) + e" 26 ' (dx + n x dy + n 2 dz) 2 + e - 2t + 2l}1 (dy + n 3 dzf , (4) 

where a ^ 1 in general. 

An important subclass of the G 2 models are the orthogonally transitive (OT) models, see [62, 63] 
and p. 43 in [7]; we refer to Appendix A.l for details. In these models the coordinates can be 
adapted so that n 2 = n 3 = and a = 1, which corresponds to setting the lapse function to 5 

N = - exp(-[t + b 3 }) = -cxpHfe 1 +b 2 + b 3 }) = -Vdet ( 3 )y . (5) 

The metric in the OT case thus simplifies considerably, see also (47) and (64). 

Throughout this paper we focus on analyzing and representing solutions by means of conformally 
Hubble-normalized Iwasawa adapted variables, see [47] and section 2 of [48], rather than the metric 
variables b , b 2 , b 3 and n\, n 2 , n^. The complete set of relations between the conformally Hubble- 
normalized variables and the metric variables are given in Appendix A. 2; here we restrict ourselves 
to giving the relations for the key objects: The Hubble scalar, which is one third of the derivative 
of the logarithm of the spatial volume density w.r.t. proper time, is given by 

H = -lN- 1 d x o(b 1 + b 2 + b 3 ), (6) 

while the conformally Hubble-normalized shear variables take the form 

T, a :=^ aa = -l-(HN)- 1 d x0 b a (a = 1,2, 3), (7a) 

£23 = \ e x P(& 3 - b 2 ) {H Nylons , £12 - \ cx P (6 2 - b 1 ) {HN^d^m , (7b) 

so that Ei + S 2 + S 3 = 0; note that £31 = 0, see Appendices A.l and A. 2. We use Greek indices 
to refer to the components of tensors w.r.t. the conformally Hubble-normalized Iwasawa frame. In 
general a spatial orthonormal frame rotates w.r.t. a gyroscopically fixed spatial frame (which is a 



5 Throughout this paper we choose the direction of time toward the past. 
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so-called Fermi-propagated frame) . In the present approach, this spatial frame rotation is encoded 
in the three functions R a , a = 1,2,3. As a consequence of the Iwasawa parametrization, these 
functions are tied to the off-diagonal shear components by 

£23 = —Ri , £12 = — i?3 , £31 = i?2 = , (8) 

see [48]. We choose to use R\ and R3 as independent variables, instead of £23 and £12. Note 
that £23 = —Ri = in the OT case, see Appendix A. 2. Finally, it is of interest to also give the 
Hubble-normalized spatial commutator function Aq (= Nn): 

Aq = H" 1 exp(6 2 + b 3 - b^dzm . (9) 

The Einstein vacuum equation in the conformal Hubble-normalized orthonormal frame approach 
are given in Appendix A. 3. 

In the following we will analyze G2 vacuum models in terms of the conformally Hubble-normalized 
orthonormal frame variables. In particular, to illustrate the behavior of solutions, we will use 
projections onto (£1, £ 2 , £s)-space and monitor the behavior of Aq along the arising trajectories. 

3 Transitions 

In the Hubble-normalized approach a state space picture emerges which is central for the analysis 
of the (asymptotic) dynamics of G2 models. The key structure on the boundary of the state space is 
the 'Kasner circle' of fixed points which represents 'generalized' Kasner metrics. These fixed points 
are connected by an intricate network of curves and families of curves representing 'generalized' 
spatially homogeneous solutions and inhomogeneous solutions, which we denote by transitions in 
this context. The Kasner circle and the transitions turn out to be the building blocks of generic 
oscillatory singularities of G2 models (and beyond); see section 5. 

3.1 The Kasner circle 

The Kasner, i.e. vacuum Bianchi type I, solution is usually given as 

ds 2 = -dP + P^dx 2 + t 2 P-df + P^dz 2 , (p! + p 2 + p 3 = 1 = p\ + p 2 + pi) , (10) 

where the constants Pi,P2,P3 are known as the Kasner exponents. By making a constant trans- 
formation of the coordinates according to 

t = c° i , x = , y — c 2 jX J , z = c 3 jX-' , (11) 

the Kasner solution takes the form 

ds 2 = -{c%dt) 2 + t 2pi (c^da?) 2 + P P2 (c 2 jdx j ) 2 + t 2p3 (c^dx 1 ) 2 , (12) 

where the constants c l j are obtained by appropriately scaling c j with c%. The generalized Kasner 
metric is defined by letting the constants p a , c°q and c l j be arbitrary functions of the spatial coor- 
dinates. As a consequence this metric is not a solution of Einstein's field equations; its importance 
lies in its role as a building block in the description of the asymptotic dynamics of actual solutions. 
In the present context of G2 models with Iwasawa parametrized metrics we have c\ = c\ = c\ = 0, 
and functions are functions of z (= x 3 ) alone. 

In Appendix B of [48] the conformal Hubble-normalized approach was used to derive the gener- 
alized Kasner metrics by perturbing away from the ' local boundary' (which previously was called 
the silent boundary, see [64, 65] for a discussion). The local boundary is an invariant boundary 
subset of the Hubble-normalized state space that is obtained by setting all Hubble-normalized 
spatial frame components (and hence spatial frame derivatives) to zero. 6 By construction, the 



In the G2 context the local boundary is obtained by setting E3 = 0, and hence E3 d z = 83 = 0, in Appendix A. 3. 
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equations induced by the field equations on the local boundary coincide with the equations de- 
scribing spatially homogeneous models, where, however, constants of integration become spatial 
functions. Accordingly, the local boundary provides a well-defined state space setting for the BKL 
concept of 'generalized spatially homogeneous metrics'. 

The structure on the local boundary that is of prime importance for our understanding of the 
asymptotic dynamics of solutions toward a singularity is the Kasner circle K°. The Kasner circle 
is obtained by setting all Hubble-normalizcd variables to zero except the diagonal shear variables 
Ei, E 2 , E3. (In particular, the rotation variables R\, R2, R3 are zero, i.e., the frame that is 
used is a Fermi-propagated frame.) The equations then imply that the diagonal shear variables 
are constants in time but functions of the spatial coordinate(s). The single remaining (algebraic) 
equation on the local boundary is the Gauss constraint, which then characterizes K°: 

E Q = const Va, E 1 + £ 2 + E 3 =0, E 2 := iE^E"' 3 = | (E? + T? 2 + £§) = 1 , (13) 

where E a are related to the generalized Kasner exponents p a according to 7 

E Q = 3p Q -lVa, pi+P2+P3 = l, p\+pl+pl = l- (14) 

The Kasner circle is divided into six equivalent sectors, which are associated with axes permutations 
and denoted by permutations of the triple (123): Sector (a/3j) is characterized by p a < pp < p-y, 
see Fig. 1. The boundaries of the sectors are six special points that are associated with locally 
rotationally symmetric (LRS) solutions, 

T Q : (E a ,E /S ,E 7 ) = (+2, -1,-1), or, equivalently, (p a ,pp,p~,) = (1,0,0) , (15a) 
Q Q : (E a ,E /8 ,E 7 ) = (-2,+l,+l) , or, equivalently, (p a ,pp,p~,) = (-§, §, §) . (15b) 

The Taub points T Q , a = 1,2,3, correspond to the flat LRS solutions — the Taub representation 
of Minkowski spacetime, while Q a , a — 1,2,3, yield three equivalent LRS solutions with non-flat 
geometry. 

A frame independent (gauge invariant) way of representing the Kasner exponents p a is by means 
of the (standard) Kasner parameter u; for each of the six equivalent sectors (otflj) we have 

p a = -u/f(u), pp = (1 + u)/f(u) , pj = u(l + u)/f(u) , (16) 

where u e (l,oo), while the boundary points of sector (a/37), Qa an d T 7 , are characterized by 
u = 1 and u = 00, respectively. The function /, which will be used extensively, see also e.g. [4], is 
defined by 

f = f(x) :=l + x + x 2 . (17) 

In addition to the standard Kasner parameter u, we follow [20] and define the extended Kasner 
parameter u by 

Pl = -u/f(u), p 2 = (l + u)/f(u), p 3 =u(l + u)/f(u), (18) 

where u e (—00,00) so that each value of u distinguishes a unique point on the Kasner circle, see 
Fig. 1. Comparing (16) and (18) we obtain a transformation between u and u for each sector: 

(123) : (1, 00) 3i = u, (132) : (0,1)3« = u' 1 , (19a) 

1 u 

(312): (-\,0)3 &=-—, (321): (-l,-i)9« = -— , (19b) 

u + 1 

(231) : (-2,-1) 3 u= — , (213) : (-00, -2) 9 u = -(u+ 1) . (19c) 



7 Constants are constants in time but functions of the spatial coordinates; in the G2 case of z. For this reason 
the Kasner exponents arc called 'generalized' exponents. 
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Figure 1: The division of the Kasner circle K° of fixed points into six equivalent sectors and six 
LRS fixed points T Q and Q Q , a = 1, 2, 3. Sector (a(3j) is defined by E Q < E^ < E 7 . The values of 
the extended Kasner parameter it, cf. (18), along K° are indicated by the values of u at the LRS 
fixed points; in addition, k = 2u + 1 is given. 

3.2 Frame transitions 

In the previous subsection the generalized Kasner metric is represented in a Fermi propagated frame 
with diagonalized shear (which is at the same time an Iwasawa frame). However, it is possible to 
represent the Kasner metric in a rotating spatial frame (e.g., a rotating Iwasawa frame), which 
leads to a time-dependent and non-diagonal Hubble-normalized shear. Since the frame invariants 
of the shear tensor, i.e., its trace, square, and determinant, must remain constants in time, the 
generalized Kasner metric in a (rotating) Iwasawa frame (for which i? 2 = 0) is characterized by 
the relations 

= Ex + E 2 + E 3 , (20a) 

1 = E 2 = i(S? + S|+El) + i( J R? + i?i) , (20b) 
const = det E Q/3 = EiE 2 E 3 - ^Rj - E 3 i?| , (20c) 

where the constant is merely a constant in time. 

A frame rotation in the 2-3-plane corresponds to R\ ^ 0, i? 3 = 0, and thus, from (20c), by means 
of (20a) and (20b), 

Ei (E 2 E 3 - Rj) = Ex (±(E 2 + E 3 ) 2 + §E? - 3) = Ei (E? - 3) = const . (21) 

Since Ex = const, this implies straight lines in (E 1; E 2 , E 3 )-space, see Fig. 2(a). A single straight 
line represents a Kasner solution (or a generalized Kasner metric corresponding to a constant value 
of u) in a frame that is rotating in the 2-3-plane. A generalized Kasner metric associated with a 
spatially dependent u corresponds to a family of lines. Analogously, for a frame rotation in the 
1-2-plane we have R\ = 0, i? 3 ^ 0, and the straight lines E 3 = const, see Fig. 2(b). We refer to the 
trajectories that individual frame rotations give rise to as single frame transitions, Tn x and Tr 3 . 
Double frame transitions Tr^, for which R\R 3 7^ 0, are not expected to play an essential role 
in the asymptotic dynamics of solutions toward a generic spacelike singularity, see [48], and we 
therefore refrain from a discussion (the interested reader is referred to Appendix B and C in [48]). 
Note that transitions (i.e., solution trajectories) are occasionally referred to as orbits. 

Linearizing the equations for R\ and i? 3 in Appendix A. 3 at K° gives 

-d R 1 \^ = (Z 2 -X 3 )\ T<p R 1 =3(p 2 -p 3 )R 1 , -d R 3 \ K3 =3( Pl -p 2 )R 3 . (22) 
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(a) T Rl (b) Tr 3 



Figure 2: Projections of the two types of single frame transitions, Trj and Tr 3 , that exist in the 
G 2 case onto (E 1; E 2 , E 3 )-space. Throughout this paper, the direction of time, as indicated by the 
arrows, is toward the past. 



Hence the frame transitions Tr x originate from sectors (321), (312), (132) and the points Q3 and 
T 2 , where R\ is unstable (w.r.t. a time directed toward the past), see Fig. 2(a). We say that a 
frame transition 7r x is 'triggered' by R\; similarly, R3 triggers the frame transitions Tr 3 in sectors 
(213), (231), (321) and at Q2, Ti, see Fig. 2(b). Frame transitions map points of K° to each other; 
these maps can be expressed as the following maps of the extended Kasner parameter u: 

T Rl : u^u" 1 , Tr 3 : (23) 



3.3 Curvature transitions 

The Bianchi type II solutions on the local boundary are characterized by the vanishing of all 
variables except for the Hubble- normalized shear variables and the commutator function N± . The 
simplest representation is w.r.t. a Fermi propagated frame in which the shear has been diagonalized 
(which is an Iwasawa frame at the same time), i.e., R a = Va. The Bianchi type II solutions can 
be conveniently parametrized as follows: 

Ei = -4 + (1 + u 2 )( , Xp=2-u 2 {, E 7 = 2-C, (24) 

where u = ?i_ G (0, +00) parametrizes the different possible initial Kasner states and £ is mono- 
tonically increasing, see (27). In the projection onto (Si, E 2 , E3)-space, (24) describes a family of 
straight lines originating from the point (Si, E 2 , S3) = (—4, 2, 2) outside of K°, see Fig. 3. At the 
same time, the variable N\ is determined by the constraint 

i (E? + E 2 + Eg) + £JV? = 1 . (25) 

We refer to the type II trajectories as single curvature transitions 7jvi • 
Linearizing the equation for N± in Appendix A. 3 at K° gives 

-do^L - -2(1 + Si)|kp Nj. = -foWi . (26) 

Accordingly we say that the curvature transitions are 'triggered' by N±, since JVi is an unstable 
variable in sectors (123) and (132) and at the point Qi (w.r.t. a past-directed time variable), see 
Fig. 3. 

The function £ in (24) is governed by the equation 



« (=8|4il I'-sJlr'J' wtere c± ^' (27) 
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Figure 3: Projections onto (Si, S2, Sa)-space of the single curvature transitions T/Vi- 



which has the solution 

C= ... U _ T 3 j7z U+T with r:=tanh(2ut), (28) 
/(u_)V" + /(«+) 2" 

so that £ is monotonically increasing from to £ + . By u = u_ we denote the initial (standard) 
Kasner parameter; u± is the initial/final extended Kasner parameter. It follows from (24) that a 
single curvature transition T/Vi gives rise to a map 

T/Vi : u = u— > m + = — it_ . (29) 

In terms of the standard Kasner parameter u this translates to u_ i— > u + = u_ — 1 if w_ > 2 and 
w_ i— > u_)_ = (u_ — 1) _1 if 1 < U- < 2, which is the usual representation of the Kasner map. 

The Bianchi type II metric assumes a less transparent form in an Iwasawa frame that is rotating 
w.r.t. a Fermi propagated frame. In this case we speak of mixed curvature/frame transitions since 
iVi 7^ and R\ or R3 7^ simultaneously, see Appendix D in [48] for a discussion. Although 
the trajectory of a mixed curvature/frame transition looks rather complicated, the final state on 
K° coincides with the point obtained by successively applying the corresponding single curvature 
transition and the relevant single frame transition(s); single transitions thus act as building blocks 
that determine the final state of multiple transitions. In any case, just like double frame transitions, 
the mixed curvature/frame transitions are expected to not play a role in asymptotic dynamics of 
solutions toward a generic spacelike singularity, see [48]. 



3.4 Spike transitions 

Together with K° the frame and curvature transitions on the local boundary constitute the essen- 
tial building blocks for the description of the asymptotic dynamics of G2 models along 'generic' 
timelines, see sections 4 and 5. Evidently, frame/curvature transitions and concatenations thereof 
concern the BKL scenario of (asymptotic) locality in which the dynamics of timelines decouple from 
their neighbors. However, to obtain a complete description of the dynamics of G2 models toward 
generic spacelike singularities, the failure of asymptotic locality along particular timelines needs 
to be taken into account. This failure of locality is associated with the formation and recurrence 
of spikes; the central building block to describe this non-local scenario is the class of non-local 
transitions ('spike transitions') associated with the 'explicit spike solutions' found by Lim [53]. 
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The explicit spike solutions are Gi OT vacuum solutions (i.e., n 2 — n 3 = 0) that in the area time 
gauge (3)-(5) take the form 8 



b 1 

b 2 
b 3 

Til 



1 

2 

t-b 1 

1 

4 
1 

2 



— t — logsech(wt) + log ^1 + (we'zsech(wt)) 2 ^ 
(3 + w 2 ) t + A log sech(wi) - 2 log (l 



we zsec 



hM)) 2 ) 



w [e~ 2t + 2 (to tanh(u;i) - 1) z 2 ] 



(30a) 
(30b) 
(30c) 

(30d) 



where w > 0. 9 



In the following we analyze the explicit spike solutions (30) in terms of Hubble- normalized variables. 
For this purpose we define 

T := tanh(ttfi) , 0(T) := (1 - T) 1_ » (1 + T) 1+ » , 5± (T,z) := 1 ± i?(T)w 2 z 2 . (31) 

The time variable T is bounded; £ e (—00,00) corresponds to T e (—1, 1). In the limit T — > ±1 
we have 

lim t?(T) = , lim 0(T) = J ° W > 1 ' (32) 



T->— 1 



and hence 



T iim_ iS± (r, z ) = i : 



T-S-+1 



T hm iff± (r,z) = i : 

W>1 



00 w < 1 . 



T hm iff± (T,z) 

U)<1 



1 z = 0, 
±00 z^0. 



Using the variable transformation given in Appendix A. 2 leads to 



Pi ■= 



P2 



P3 



+ 

(1 + S 2 ) 
(I + S3) 



2 W (T (7» 



(w 2 + 3) ff4 


(T,z)- 


Aw (T- 


--) ' 


4.9+ (T, z)- 


- 2w (T 


w / -J 


(T,z) 


(w 2 + 3)g. 


f(?>)- 


-Aw (T 


-i) 

W J 


( w 2 - 1) 3+ 


(T,z)- 


Aw (T - 




( W 2 + 3)ff4 


{T,z)~ 


Aw (T - 


--) ' 



where pi + P2 + P3 = 1, since Si + S 2 + S3 = 0, and where 



N 1 = 12w tt^tPi 

9-{T,z) 



Vl -T 2 
^3 = -3 r _ y pi ■ 



(33) 

(34a) 
(34b) 
(34c) 

(35) 



From these expressions it is straightforward to describe the behavior of the explicit spike solu- 
tions (30). Henceforth, by 'explicit spike solution' we typically mean its representation in Hubblc- 
normalized variables and its projection onto (Si, S 2 , S 3 )-space. 

In contrast to frame and curvature transitions, an explicit spike solution consists of an entire 
family of curves parametrized by the spatial coordinate z. As t — > —00 (T — »■ —1), each spike 
solution (30) converges to a point on the Kasner circle K° described by the extended Kasner 
parameter u = u_ e (0, +00), which we call the initial extended Kasner parameter, given by 



u = it- = - (w + 1) 



(•» w = 2u - 1) 



(36) 



8 The solutions are presented here in a slightly different form than in [53, 54]. Note also that, as stated in [54], 
there is a typographical error in equation (34) in [53]. For further discussion, see section 4. 

9 It is possible to let w £ R; however, w 1— > — w has the same effect as the coordinate reflection y 1— > — y and 
thus results in the same solution, cf. (42); in addition, the solution with w = is trivial and not associated with a 
transition since it corresponds to a Kasner fixed point with u = 1/2. 
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This relation becomes manifest when we use an alternative form of (34) replacing w by 2u — 1, 

g _ -[«(^) + (i-«)m] g -(r,z) 

Pl /(n)(i^)+/(l-ti)(i±^)+/(n-l)i?(TK 2- ^ i 

Setting T = — 1 and comparing with (18) establishes u as the initial extended Kasner parameter 
according to (36). Note that the requirement w > corresponds to u > |. If u > 1, it follows that 
the initial sector is (123), while if | < u < 1, the initial Kasner solution is in the part of sector 
(132) with 1 < u < 2 (where we recall from (19) that u = u if u > 1 and u = if < u < 1). 

In the limit t — > oo (T — > 1), however, a qualitative difference arises that leads to a classification 
of spike solutions in the OT G2 context: 

Hi A spike solution with w > 1 (i.e., u — U- > 1) yields a Tug/i velocity' spike transition, which 
we denote by Th\- (The term velocity is defined in Appendix A. 4.) Its limit as t — > —00 
(T — > —1) is a Kasner point in sector (123) with u = u_ = (w + l)/2, while its limit as 
t — > 00 (T — > 1) is a Kasner point that belongs to one of the four sectors (2(3j), (3(3j). This 
final Kasner point for the spike solution is characterized by an extended Kasner parameter 
u + . Letting T — > 1 in (34a') shows that 

T Hi : u+ = 1 - u- . (37) 

By means of (19), the map (37) can be expressed in terms of the standard frame-invariant 
Kasner parameter according to the map: 




e [3,oo), 



Tur. u + = l)r - ^ U - S g' 3 J' (37') 



- l)- 1 - 1) [3/2,2] 
«- G [1,3/2], 

where u — u_ and u + are the initial and final Kasner parameters, respectively. We refer to 
Figs. 4(a) and 4(b) for a depiction of the projection of high velocity spike transitions onto 
(Si,S 2 ,S 3 )-space. 

Lo A spike solution with < w < 1 (i.e., | < u = u_ < 1) yields a 'low velocity' spike solution, 
which we denote by S\_ D . Its limit as t — > —00 (T -4- —1) is a Kasner point in the part of 
sector (132) with u > \ (u < 2), which yields u = u_ = 2/(w + 1) since u = vT 1 . Its limit as 
t — >• 00 (T — > 1) is a generalized Kasner metric with discontinuous Kasner parameter which 
is because the limit of g±(T,z) as T — > 1 is discontinuous, cf. (33). From (34a') it follows 
that 

fl-U_ 21 = 0, 

5 Lo : fi+(«) = r , ^ n (38) 
^it_ — 1 z 7^ , 

which in terms of the standard Kasner parameter results in 



(u-(u--l)- 1 z = 0, 
u + {z) = \ { (38') 
I (u_ — 1) z 7^ . 



Therefore the limit (as f — > 00) of the low velocity spike solution along the timelines of 
the spike surface z — differs from the limit along the other timelines. The trajectory in 
(Ei, E2, Es)-space representing the spike surface ends in the part of sector (132) with u > 2 
(0 < u < 5), while trajectories associated with timelines with z ^ meet in a point of sector 
(312); see Figs. 4(d) and 5. This behavior of the low velocity spike solution results in the 
convergence of, e.g., the Hubblc-normalized Kretschmann scalar to a discontinuous limit. 

The trajectories of an explicit spike solution that represent timelines with \z\ ~> 1 (i.e., timelines 
far from the spike surface z = 0) are approximated by a succession of (two or three) curva- 
ture/frame transitions, see Fig. 4. (In sections 4 and 5 we will refer to finite and infinite sequences 
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of curvature/frame transitions as BKL chains.) In the case of a high velocity spike transition Th\, 
trajectories with \z\ ^> 1 are approximated by the chain Tzvi — Tr 3 — Tni\ in the case of a low 
velocity spike solution <S|_ OJ \z\ ^> 1 trajectories are approximated by a Tn x followed by a Tr 3 tran- 
sition. Formally, these chains are the limit as z — »■ ±oo of the orbits representing a spike solution. 
Since T/Vi is associated with the map u i-> — u while Tr 3 gives it >->• — (u + 1), it follows that 

Hi: u_ 2^ ^ + 1) 2^ 1 - u_ = u + (39) 

in the high velocity case, which coincides with (37) . In the low velocity case the final T/Vi transition 
is missing so that 

Lo : u_ -u_ -(-u_ + 1) = u_ - 1 = ii + , (40) 
which coincides with the z ^ case of (38). 

The evolution of the Hubblc-normalized variables along the timelines of the spike surface z = 
itself is represented by a single straight line in (£i, E2, E3)-space, which we refer to as the 1 spike 
surface trajectory' Ts- The family of 'spike surface trajectories' (parametrized by w or u_) has a 
common focal point: Each spike surface trajectory Ts corresponds to a straight line that originates 
from the point (E 1; £ 2 , £3) = (— |, \, 2) outside K°, see Fig. 6. In the high velocity case, the spike 
surface trajectories enter the physical state space in sector (123), while in the low velocity case, 
they enter the physical state space in the part of sector (132) with 1 < u < 2, i.e., \ < u < 1. 
From (34) we see that 

SiU=o =-§ + |/H*±)x, £ 2 Uo = l+X- \f{-u±)x, £ 3 | Z=0 =2- X , (41a) 

where 

X = /( fi _)^+/(fi + )l^ ' (41b) 

so that x increases monotonically from X- = 3//(u_) to x+ = 3/f(u + ). The expression /(— tt±) 
appearing in (41a) stands for /(— w_) or /(— either can be used since they coincide. 

Finally, the trajectories representing timelines z 7^ that are not far from the spike surface 
'interpolate' between the two types of behavior, i.e., \z\ » 1 and z = 0, described previously. 
This family of trajectories is rather complicated, see Fig. 4 and [53]; we note, however, that all 
trajectories intersect at T = 1/w in (Ei, £2, Ss)-space, as follows from (34). 

In equation (30) we have presented the principal form of the explicit spike solution. There exist, 
however, entire families of explicit spike solutions with different coordinate widths and time offsets. 
These spike solutions are required when spikes are concatenated, which we discuss in sections 4 
and 5. We therefore consider the constant coordinate transformation (x, y, z) 1— > (x, y, z) given by 
x = c\ x + c\ y, y = c 2 2 y, and z = c% z; in addition we may consider a rescaling of time, i.e., 
t t with t = c° i. This coordinate transformation preserves the form (47) of the metric provided 
that c°o = |c 1 iC 2 2C 3 3| and induces the transformation (6 1 , b 2 , 6 3 , m) \—¥ (b 1 , b 2 , b 3 , Hi) with 

b 1 ^b 1 -\og\c\\, b 2 = b 2 - \og\c 2 2 \ , b 3 = b 3 -log\c\\, ni = ^ + ^m. (42) 

Choosing c\ = (2Q y 1/2 , c\ = (2Q )- 1/2 Q 2 , c\ - (2Q ) 1/2 , and c 3 3 = exp(-A 2 /4) reproduces 
the free constants Qo, Q2, and A 2 of Eqs. (33)-(35) in [53] and Eqs. (A3)-(A5) in [54]; note, 
however, that {c\, c\, c\, c%} are four free constants. Restricting ourselves to the case c°o = 1 
(so that [c^c^c^l = 1) it follows from conceptual considerations or directly from the relations of 
this section and Appendix A. 2 that the Hubble scalar and the Hubble-normalizcd variables are 
unaffected by (42), although the frame variable E 3 , cf. (58), transforms like E 3 1— > E 3 = (c 3 3 ) _1 E 3 . 
As a consequence, equations (34) and (35) take the same form, where z i-> c z z, c z := c%, so that, 

6 § 

g ± (T,z) = l±i>(T)w 2 c 2 z z 2 . (43) 

Equation (30) with z ^ c z z yields a family of spike solutions and transitions with different co- 
ordinate widths, an inverse measure of which is, e.g., the constant c z . It follows from (35) that 
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(a) Tm with u- = 2.25 (w = 3.5) (b) T H ; with u_ = 1.25 (w = 2.5) 




(c) Borderline case: u— = 1 (w = 1) (d) S|_o with ii_ = 0.75 (w = 0.5) 



Figure 4: High velocity spike transitions Th\ and low velocity spike solutions <S|_o in their projections 
onto (Ei, E2, E3)-space. Each Th\/S\_o corresponds to a family of curves parametrized by the spatial 
coordinate z. The thick (straight) lines represent the trajectories of the timelines of the spike surface 
z = 0, i.e., spike surface trajectories Ts- The thin curves represent the trajectories of timelines 
7^ \z\ 1. The thin (straight) lines represent \z\ ^> 1 timelines; these are short BKL chains, 
i.e., sequences of curvature/frame transitions. The dashed lines represent constant time slices of a 
solution. 



c z = (l2itfVTOpi|«=o) d z N x | 2=Q . (44a) 

To eliminate T we make use of the trajectory of the spike surface z = 0, because this trajectory is 
not affected by the spike width; using pi\ z =o, see (34a), we obtain 

HT) = 2n t4 1 H+ - /("+)PiU=o) 1_ - (u- +/(fi-)pi|,=o) 1+ " , (44b) 

where u_ (w + ) is the initial (final) point of 7hii or 5l (for z = 0), i.e., = 1 — u_, w = 2u_ — 1. 
Inserting (44b) into (44a) we obtain an explicit formula for c z in terms of certain Hubble-normalized 
variables (and a spatial derivative). 

While a rescaling of the spatial coordinate z yields a spike solution with a different coordinate width, 
a translation in z puts the spike surface at a different spatial location. Analogously, translations 
in time ('time offsets'), i.e., 1 1— > t±to, yield spike solutions with a different temporal flow (which, 
of course, does not affect the orbits in the diagrams). In this case, the 'time' T of (31) is redefined 
as T = tanh[w(t ± t )]. 
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Figure 5: The 'skeleton' of a <S|_o spike solution described by the spike surface trajectory 7s (thick 
line) and the short BKL chain representing the trajectories associated with \z\ 3> 1; cf. Fig. 4; in 
addition the dashed line Sa corresponds to a 'missing' curvature transition Tn-l that would reunite 
the spike surface z — with the z ^ trajectories. 




Figure 6: Projections of spike surface trajectories Ts onto (Si, S 2 , S 3 )-space. 



3.5 Bill spiky features 

High velocity spike transitions %\\ originate from Kasner points with 1 < u < oo, i.e., from sector 
(123); low velocity spike solutions Si D from points with | < u < 1, which covers merely the upper 
part of sector (132). The third type of spiky feature, which we call a Bill spiky feature, fills the 
rest of sector (132): 



Bill The limit of a Bill spiky feature as t — > — oo is a Kasner point on sector (132) with < u < 



l 

2 ' 

i.e., u > 2. The trajectories z ^ of a Bill spiky feature coincide with the 7tvi transition 
that originate from the Kasner point u on sector (132) and end in sector (312), while the 
spatial points of the spike surface z — are left behind at the Kasner point u on sector (132). 
Accordingly, the limit as t — > oo is a generalized Kasner metric with discontinuous Kasner 
parameter such that 

fii- z = , 
— u_ z f= . 

In terms of the standard Kasner parameter we have u_ h- > u + = u_ and u + = it_ — 1, 
respectively. 
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A Bill spiky feature is intimately connected with the 'second half of a low velocity spike solution 
<S|_ m the neighborhood of z = 0. Specifically, consider a S\_ characterized by w — 1 — 2u with 
u G (0, 5), (inverse) coordinate spike width c z , and a time offset t , i.e., t is replaced by i + t 
in (30). For a given value of z = z 7^ 0, if c z is sufficiently small, the trajectory of z will resemble 
the concatenation of the low velocity spike surface trajectory 7s and a T/Vi transition originating 
from u. It is possible to choose a (large) value of to, which is fine-tuned w.r.t. c z , such that the 
point t = of the trajectory is close to the midpoint of the Tn 1 transition. This motivates taking 
the combined limit c z — > and to — > 00 such that 2c z e( 1 ~ u ')* = 1, which, via (34), yields 

1. , - (1 - 2t2) 2 z 2 e 4 "*) , 

Pl = 3 (1 + Sl) = /(«)+/(-«)(! -2fi)^ e «t ■ ( 46a ) 



1 , ^ , _ 2 - (1 - u) (1 - (1 - 2u) 2 z 2 e 4 " 4 ) 
/(i2)+/H2)(l- 2*2)2 



p 2 = -(1 + S 2 ) = ^ , ^ ^ o.-.^^t . ( 46b ) 



-12u(l - 2u)ze 2 "' 

1 ~ f(u) + f(-u)(l-2u) 2 z 2 e 4M ' ( ' 

while i? 3 and Aq 2 are zero. The above expression is not a solution to Einstein's equations, but it 
is an inhomogeneous solution of the equations on the local boundary. 



4 Asymptotics and concatenation in the OT case 

In the present section we restrict our attention to the OT subclass of the G2 models, of which 
the T 3 Gowdy models are a prominent special case, see Appendix A.l. The metric of the OT G 2 
models is of the form (1), where n 2 = n 3 = 0, i.e., 

ds 2 = -N 2 (dx°) 2 + e- 2bl (dx + ni dy) 2 + dy 2 + e" 26 '^ 2 , (47) 

where b a = b a (x°, z) \fa and n\ = ni(x°, z). The area time gauge entails x° = t = b 1 + b 2 (modulo 
a constant) and N = —y/dct (%. Recall from section 2 that the conformally Hubble-normalized 
rotation variable R\ vanishes identically under the OT assumption. 



4.1 BKL concatenation in OT compatible Bianchi models 

Let us briefly recapitulate the basic facts about spatially homogeneous models in the present 
context. The existence of two commuting Killing vector fields is incompatible with the defining 
properties of the Bianchi type VIII and IX models. The Bianchi types that are compatible are 
types VII/,, and VI^, and the 'lower types', which can be obtained by means of Lie contractions. 
However, it is only the symmetry adapted frames of Bianchi types VIIo, VLj, II, and I that are 
compatible with a line element (47) in the area time gauge, cf. [61]. 

The initial singularity of the spatially homogeneous OT G 2 vacuum models is known to be a 
Kasner type singularity, i.e., the asymptotic behavior is that of a Kasner solution. This can be 
understood in connection with Fig. 7. For Bianchi models of type VI and VII , the non-trivial 
Hubble-normalized variables are the diagonal shear variables and N\ and R3. From (22) and 
Fig. 2(b), and (26) and Fig. 3, respectively, we see that R3 is an unstable variable ('triggered') in 
three sectors of the Kasner circle while Ni is triggered in two sectors, see Fig. 7(a). Sector (312), 
however, is a stable sector, which means that Aq and R3 are stable variables, i.e., a solution with 
initial data close to a Kasner point of sector (312) must converge to a Kasner point of that sector, 
because Aq and R 3 decay exponentially. It is thus intuitive (and straightforward to prove, see, 
e.g., [48]) that every type VIo or type VIIo model converges to a fixed point on the stable sector. 

Considering Figs. 2(b) and 3 together, it is clear that frame transitions and curvature transitions 
can be concatenated by identifying the 'final' fixed point (cj-limit point) of one transition with 
the 'initial' fixed point (a- limit point) of another transition, see Fig. 7(b). Concatenation of Tr 3 
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and T/Vi transitions yields finite heteroclinic chains, 10 which we refer to as OT BKL chains; these 
chains terminate on the stable sector (312). 




(a) OT triggers (b) OT BKL chain 



Figure 7: In the OT case each sector is associated with a 'trigger' {N\ or i? 3 ) except for sector 
(312) at which these variables are stable. As a consequence, OT BKL chains (i.e., concatenations 
of T/Vi and Tr 3 transitions) are finite and terminate at sector (312). 

The relevance of the OT BKL chains is immediate: The evolution of initial data sufficiently close 
to a Kasner point of any other than the stable sector sector (or: the evolution of data sufficiently 
close to a transition) is approximately described an OT BKL chain, which yields a finite number 
of 'oscillations' between Kasner states, mediated by Tni an d Tr 3 transitions, and convergence to 
a Kasner point on the stable sector (312). However, it should be pointed out that the evolution 
of initial data far from the Bianchi type I and the type II subsets (i.e., far from the Kasner circle 
and the transition orbits) need not bear any relation to OT BKL chains, although there is still 
convergence to a Kasner point on the stable sector. For a comprehensive study of type VIo/VIIo 
dynamics we refer to [66, 41] and references therein. 

4.2 Asymptotics of spatially inhomogeneous OT models 

According to the BKL conjecture, asymptotic locality means that in the inhomogeneous case there 
exists a generic set of timelines whose evolution is governed asymptotically by the flow on the local 
boundary (which entails that the evolution of a single timeline is congruent with the evolution of 
a spatially homogeneous model.) This does not mean, however, that the dynamics of spatially 
homogeneous models is sufficient to describe the dynamics of inhomogeneous models. On the 
contrary, spatially inhomogeneous OT models behave quite delicately in general. 

However, Fig. 7(a) suggests a straightforward behavior for at least certain classes of initial data. 
Consider, e.g., initial data that constitute a small perturbation of a generalized Kasner metric 
associated with the stable sector (312) on K°, i.e., Si, £2, S3 are functions of the spatial variable 
z with values close to sector (312), while the remaining variables and the Hubble- normalized spatial 
frame derivatives are small; note that we do not exclude that, e.g., Ni has zeros. Since neither 
Ni nor i?3 are triggered but decay exponentially, we expect convergence to a generalized Kasner 
solution on sector (312). This is an asymptotic state that is often termed 'asymptotically velocity 
dominated'; we refer to Appendix A. 4, where we discuss results that turn the intuitive picture into 
a rigorous one in the case of T 3 Gowdy models. 

Prescribing 'almost Kasner' initial data close to sector (321), on the other hand, leads to an 
exponential decay of N±, while R3 is triggered, so that an (approximate) Tr 3 frame transition 



10 A heteroclinic orbit is an orbit (i.e. solution trajectory) that starts («-limit) and ends (t»;-limit) at two different 
fixed points. A heteroclinic chain is a sequence of heteroclinic orbits such that the w-limit point of one orbit is the 
Q-limit point of the subsequent orbit. 
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ensues at each spatial point, 11 which in turn leads to the convergence of the solution to a generalized 
Kasner metric of sector (312), as before. 

Likewise, initial data close to sectors (231) or (213) lead to an initial frame transition 7~r 3 , which 
takes the solution to sectors (132) or (123) where N± is the trigger. If |A^| > (in the optimal case: 
uniformly) we are led to expect the familiar behavior of the spatially homogeneous type VIo / VII 
case, i.e., a subsequent (approximate) T/Vi curvature transition at each spatial point; accordingly, 
each timelike follows an (approximate) OT BKL chain. On the other hand, if N\ has zeros for 
certain values of the spatial coordinate z, the simple picture inspired by spatially homogeneous 
dynamics fails. At these spatial surfaces it cannot be the growth of Ni that drives the further 
evolution of the solution. Instead, spatial gradients take the role of triggers at these special values 
of z (and their neighborhoods) and spatial structures ('spikes') develop (which are 'true' spikes 
in the sense of being gauge invariant features). These spikes are approximately described by the 
explicit spike solutions 7hi, Slo, and the Bill spiky feature in a sense that we discuss in section 4.3. 

We emphasize that the evolution of general initial data is rather complicated (and the heuristic 
considerations fail). However, asymptotic velocity dominance remains true in general, i.e., solutions 
converge to a generalized Kasner metric (which is in general associated with a discontinuous Kasner 
parameter); see Appendix A. 4. 

4.3 Spike concatenation in OT models 

In section 4.1 we have concatenated curvature transitions and frame transitions to build OT BKL 
chains. Analogously, high velocity spike transitions Th\ and frame transitions Tr 3 can be concate- 
nated by identifying the 'final' Kasner point of one transition with the 'initial' Kasner point of 
another transition, see Figs. 8(a) and 8(b). Note that this is possible because the entire family of 
curves representing a Th\ transition converges to a point on K° as t — »■ — oo and to another point 
on K° as t — > oo. Concatenation yields finite 12 chains, which we refer to as high velocity chains. 

Tracking the extended Kasner parameter along a high velocity chain we find an alternation between 
the 7hi law u n- 1 — u and the Tr 3 law ii4 — (u + 1), whose combination results in u n- u — 2. 
A high velocity chain either ends at a point on the stable Kasner sector (312), where the final 
transition can be a 7r 3 as in Fig. 8(a) or a 7h\ as in Fig. 8(b), or the chain ends at a point on 
sector (132). The latter case is particularly interesting, since the high velocity chain can then be 
continued by a <S|_o solution or a Bill spiky feature, depending on its location in sector (132). We 
use the term 'OT spike chains' as an umbrella term for the two cases, i.e., an OT spike chain 
is either a high velocity chain ending in sector (312) or a high velocity chain followed by a <S|_ 
solution or a Bill spiky feature. Note that the continuation of OT spike chains beyond this point is 
impossible (because Tr x transitions do not exist in the OT context; see, however, the next section). 
Note that for an OT spike chain ending with a S\_ a solution or a Bill spiky feature the trajectories 
of the points on the spike surface z = end at sector (132) with < u < 1/2, while points with 
z^O are transported to the stable sector (312). Therefore, in a vague sense, the part of sector 
(132) with < u < 1/2 is 'stable' for the spike surface, see Fig. 9. 

The role of OT BKL chains for Bianchi models of type VIo and VIIo has been discussed in sec- 
tion 4.1. Numerical evidence [54] suggests that OT spike chains fill the same function for inho- 
mogeneous solutions with spikes. In the following we give a detailed description. For simplicity, 
we restrict our attention to symmetric initial data sets, by which we mean initial data such that 
Si, £ 2 , S 3 , and R 3 are even functions in z, while Ni and N X2 are odd functions (where we are 
interested in the non-degenerate case for which the gradients at z — are non- vanishing). 13 This 

11 Technical difficulties occur if R3 is zero initially at some value of z; in this case a 'false spike' ensues. We 
refrain from discussing this case since it is a gauge effect (a "nuisance" as stated in [67]) and focus on 'true' spatial 
structures instead. We note however that a complete description of a generic singularity in an Iwasawa frame 
necessarily involves gauge features of this kind. 

12 Concatenation in the opposite direction of the time coordinate, i.e., away from the singularity, yields infinite 
chains converging to the Taub point T3. 

13 In addition we assume the 'genericity condition' that R3 ^ at the spike surface 2 = 0. 
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(c) Th-Tr 3 -S Lo (d) Th -T R3 -Th-Tr 3 -B\\\ 



Figure 8: OT spike chains in the projection onto (Ei, E2, £3)-space; depicted are chains of the 
spike surface trajectories (z = 0) and trajectories associated with \z\ ^> 1 (which are OT BKL 
chains). An OT spike chain is a high velocity chain (i.e., a concatenation of 7~h\ transitions and 
7~r 3 frame transitions) ending with a Tr 3 transition (a) or a 7~h\ transition (b); or there is one 
additional (final) transition: a low velocity solution S\_ (c) or a Bill spiky feature (d). The latter 
is depicted by a dashed line. 

type of initial data fixes the location of the spike surface at z = 0. 14 Note that the explicit spike 
solutions Tw, S\_ , and the Bill spiky feature, see (34) and (46), respectively, are of this symmetric 
type. Since the location of the spike surface is fixed, it is simple to calculate the particle horizons 
associated with the spike world sheet z — 0. Since a null vector orthogonal to the surfaces of 
symmetry is of the form (1, 0, 0, ±e~') T , see (47), the particle horizon at t of the z — timelines is 
[— e _t ,e _t ] times the symmetry surfaces; in particular, the particle horizons are shrinking rapidly 
(in z-direction) as t — > 00. Hence the dynamics of a solution in a neighborhood of the spike world 
sheet for t > to is completely determined by the initial data of the solution in a neighborhood 
of [— e~ to , e _t °]. In particular, changing the initial data for an explicit spike solution outside a 
neighborhood of [— e~ to ,e -t °] does not affect the evolution (30) of the spike surface and its im- 
mediate neighbors. 15 In the following we restrict our attention to (symmetric) initial data on a 
neighborhood of [— e~*° , e~*°] and consider the evolution of that data on its domain of dependence. 

A solution S(t,z) with initial data sufficiently close (in, e.g., a C k,a norm) to an explicit spike 

14 In general, this type of initial data fixes the location of a spike surface at 2 = 0. During the evolution of the 
data, it is possible that other (in general 'moving') spikes form at different values of z. We refer to the concluding 
remarks. 

15 The asymptotic dynamics is thus local in the sense that global topological aspects are irrelevant once the particle 
horizon scale has become sufficiently small compared to the 'global scale' associated with the topology. 
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Figure 9: 'Triggers' for the evolution of the spike surface z = in the OT case. In sectors (213), 
(231), and (321), i? 3 triggers frame transitions Tr 3 . Since Ni — at z = 0, Ni does not exist 
as a trigger; instead spatial gradients ('£') take the role as triggers and lead to a spike surface 
trajectory 7s, see Fig. 6. These triggers cover sector (312) and (132) with u > 2. Sector (312) and 
sector (132) with u < 2 are stable in this context. 

transition Tu\ will remain close to that solution for some time interval; temporarily we may thus 
speak of an approximate spike transition. The solution Sit, z) is approximated by 7~h\ at least up to 
a point when the two solutions reach a neighborhood of the Kasner circle K° . While T\\\ continues 
to converge to a fixed point on K°, the solution S(t, z) will, in general, stray off course — it ceases 
to be an approximate spike transition. The initially small deviation of S(t, z) from Th\ amplifies (in 
the sense of increasing relative errors) and leads to different behavior in the variable R3. In contrast 
to 7(Hi, this variable is triggered at K° and S(t, z) undergoes an approximate Tr 3 frame transition. 
This behavior can be understood, qualitatively and quantitatively, by considering a linearization 
of the equations at K°, but since such an analysis goes beyond the scope of the present paper it 
will be pursued elsewhere. For the present purposes we restrict ourselves to referring to the strong 
numerical evidence we have obtained for the described behavior: The numerics suggest that the 
solution S(t, z) shadows an OT spike chain 



and thus S(t, z) is an approximate OT spike chain. Furthermore, the numerical evidence shows 
that an approximate OT spike chain S shadows an exact OT spike chain with an increasing degree 
of accuracy; in particular, the asymptotic limit of S(t, z) is a small perturbation of the limit of the 
OT spike chain itself. 

For an approximate spike transition we define its (inverse) spike width c z through (44), i.e., 



will in general not be a constant exactly, but depend on when (49) is evaluated. However, for an 
approximate spike transition, c z will be approximately constant, at least up to the point where 
the deviation from Tu\ becomes too large. Note that definition (49) is not unique in the sense 
that it can be replaced by similar definitions involving, e.g., N12 and pi instead of N\ and p\. 
Evidently, for exact spike transitions this does not cause any difference, while for approximate 
spike transitions the resulting value will indeed differ from (49) in general. However, the slight 
ambiguities quickly converge to zero when the approximate spike transitions are approximated by 
the exact transitions with an increasing degree of accuracy. 




Bill 



(48) 
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Considering an approximate OT spike chain S(t, z) a natural question to ask concerns the behavior 
of the evolution of the spike width. Suppose that S(t, z) is approximated for some time interval by 
an exact spike transition 7h\ with parameter Uq, inverse width cq, and time offset — to, i.e., by (30) 
with the replacements z i— > cqz and t t — t$. In a (much) later time interval, the solution S(t, z) 
is then approximated by the subsequent Th\ of the chain; we have u\ = uq — 2, inverse coordinate 
width ci, and time offset t\ (where t\ ^> to). In Appendix A. 5 wc show that 

c x =e tl -* o c = e At c . (50) 

The time At between two subsequent (approximate) spike transitions is dominated by the time 
the solution spends in the neighborhood of K°. The closer the solution S(t, z) is to an exact OT 
spike chain, the greater At will be. As a consequence the coordinate width of the spike decreases 
rapidly toward the singularity, since the inverse width increases rapidly. 



5 Concatenation in the general G2 case 

In contrast to the OT case, the Hubble-normalized frame variable R\ does not vanish identically 
for the general Gi models. As a consequence, sector (312), which is stable for the OT models, see 
Fig. 7(a), possesses an unstable manifold which is associated with the variable R\, see Fig. 10. The 
same is true for the part of sector (132) with < u < 1/2, which is 'stable' for the spike surface 
trajectories in the OT context, see Fig. 9, but unstable through Tr x frame transitions for general 
G2 models, see Fig. 10(b). Intuitively it follows that the instability of the entire Kasner circle 
(except for, possibly, the Taub points) prevents convergence of (generic) solutions to generalized 
Kasner solutions and leads to oscillatory singularities. In this section we present analytical and 
numerical evidence that corroborate this conjecture. 




(a) Triggers for BKL chains (b) Triggers for spike chains 

Figure 10: Triggers for general G2 models. Subhgure (a) and (b) generalizes Fig. 7(a) and Fig. 9, 
respectively. 



5.1 BKL concatenation and BKL chains 

Concatenation of frame transitions Tr x , Tr 3 , and curvature transitions 7/Vi yields heteroclinic 
chains, which we refer to as a l BKL chains', see Fig. 11. As opposed to OT chains, BKL chains 
are infinite in general. 

We note that a BKL chain is not uniquely determined by an (initial) fixed point on K°. 16 This is 
because there are two sectors with two triggers (i.e., unstable modes) instead of one, see Fig. 10: 

16 Note that this is different from the Mixmaster case associated with Bianchi types VIII and IX. Recall that the 
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Figure 11: A BKL chain is an (in general) infinite heteroclinic chain consisting of T/Vd and Tr 3 
curvature and frame transitions. The paths are not unique; e.g., at the point X two continuations 
are possible, a frame transition Tr x or a curvature transition 7~/Vi (dashed line). However, the 
different paths rejoin at the point Y after an additional curvature or frame transition. At the 
point Y the chain continues with a (not shown) Tr 3 transition. 

In sector (132) it is possible to continue the chain with a Tn ± curvature transition or a Tr 1 frame 
transition; in sector (321) it is possible to continue with a Tr x or a Tr 3 frame transition. Note, 
however, that different paths rejoin after two or three transitions, see Fig. 11; in any case, the 
scenario of two triggers at a fixed point requires closer attention; we refer to the discussion below. 

In analogy with the role of OT BKL chains in the OT context, the relevance of BKL chains in 
the present context is twofold. On the one hand, the role of BKL chains for spatial homogeneous 
models is immediate: Since the symmetry-adapted frames of general Bianchi type VI_i/ 9 models 
are Iwasawa frames, these models are of the form (1) where the metric is manifestly spatially 
homogeneous, see [48]. The asymptotic evolution toward the initial singularity of general type 
VI_i/g models is expected to be represented by 17 BKL chains which lead to infinite oscillations 
between different Kasner states; the singularity can therefore be referred to as being oscillatory. 
On the other hand, there exists heuristic evidence [48] and numerical support [50, 52] that the 
asymptotic evolution towards a generic spacelike singularity along a generic (open) set of timelines 
of a general inhomogeneous vacuum (or 'vacuum dominated') solution of the Einstein equations 
(with \NiRiR^\ > 0) is represented by BKL chains as well. (For the perfect fluid case and 'vacuum 
dominance' we refer to [65] and references therein.) 

Asymptotic shadowing necessarily brings solutions (i.e., type VI_i/g solutions or single timelines of 
inhomogeneous solutions) to a neighborhood of a Kasner point where two variables are 'triggered' 
(i.e., grow exponentially). Consider, e.g., the case of a trajectory following a Tr 3 transition to 
a fixed point X on sector (132) of K°, where Ni and i?i are triggers, see Figs. 10 and 11. If 
\Ni\ <C I-Rij (<C 1) when the trajectory enters a neighborhood of the fixed point X, R\ will be 
triggered and a Tr 1 will follow; if \R\\ <C |iVi|, N\ will be triggered and a T/Vi will follow. Clearly, 
it is possible to prescribe initial data such that both R\ and N\ grow to macroscopic size, which 
means that a so-called multiple transition (here: a mixed curvature/frame transition) follows. 

symmetry adapted frames of types VIII and IX are not of the Iwasawa type; instead it is strongly preferable to use a 
symmetry group compatible Fermi propagated orthonormal frame, which yields three types of curvature transitions 
associated with the three spatial directions, while frame transitions are absent. The associated BKL chains do not 
exhibit the ambiguities of Fig. 11; see, e.g., [42] for details. 

17 The exact meaning of 'being represented by' is a delicate point. Clearly, the evolution of initial data close 
to a Kasner point is approximately described, at least for a finite time, by a finite piece of a BKL chain ('finite 
shadowing'). However, the study of Mixmaster dynamics [43, 44, 45] suggests the possibility that a generic solution 
asymptotically follows 'its' BKL chain forever. 
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Generic BKL chains enter a neighborhood of the fixed point X infinitely many times, and generic 
solution trajectories are thus expected to do likewise. However, the analysis of [48] suggests that 
the multiple transition scenario does not occur infinitely often; generically, in the asymptotic limit, 
trajectories follow single transitions (and thus BKL chains), a result that is consistent with the 
billiard picture of Damour et al [49]. The question of which of the two possible transitions, Tr x or 
T/Vi ! is followed more frequently, is delicate. Linearization the system at X yields the eigenvalues 
associated with the variables R\ and Ni . If the eigenvalue associated with, say, Ri is considerably 
larger than the other, triggering a Tr x transition seems more likely than a T/Vi transition. It is 
possible, however, that R\ is almost always sufficiently small against N\ when the trajectory comes 
close to X, which would make T/Vi transitions more frequent than expected. In any case, the fact 
that different paths rejoin, makes the ambiguity in paths a transient problem. 

5.2 Spike concatenation in the general G 2 models 

While the oscillatory asymptotic dynamics of generic timelines is represented by BKL chains, the 
evolution of spike timelines, and their neighbors, bears no relation to spatially homogeneous dy- 
namics; the BKL picture breaks down. This is because N\ is zero at the spike surface which results 
in a loss of the trigger iVi in sectors (123) and (132); in other words, Tn-i curvature transitions 
cannot take place. However, from Fig. 10(b) we see that the general G2 case has at least one trigger 
at each sector of K° (which is due to the existence of a non-zero i?i), where spatial gradients take 
the role as triggers in sector (123) and part of sector (132). We will argue that this 'instability' of 
the Kasner circle results in oscillatory dynamics which is represented by infinite spike chains. 

As in the case of OT spike chains, see section 4.3, high velocity spike transitions Th\ and Tr % & Tr 3 
frame transitions can be concatenated straightforwardly since the entire family of curves re- 
presenting a TPii converges to a point on K° as t — > —00 and to another point on K° as t — > 00. 
However, the situation for low velocity spike solutions <S|_o is different, because the limit on K° of 
the trajectory of the spike surface z = differs from the (common) limit of the z ^ timelines, 
see Figs. 4(d) and 5. It is obvious that low velocity spike solutions S\_ a do not immediately fit into 
the network of transitions. 

However, we will see that low velocity spike solutions naturally combine with a frame rotation in 
Ri and part of a high velocity transition to form a joint low/high velocity spike transition which 
we denote by 7j - This spike transition is characterized by the same fundamental property as Tni 
transitions that the entire family of curves representing 7j emerges from one point and ends at 
another point on K°. There is thus a stark contrast between the OT case and the general Gi case. 
In the OT case, a S\_ solution and perturbations thereof develop 'permanent' spatial structures 
since the timelines of the spike surface and the non-spike timelines converge to two different Kasner 
points, see Figs. 4(d) and 5; in the general Gi case, these spatial structures form and 'un-form' 
since S\_ a is merely the first part of the transition 7j : The timelines rejoin eventually. Spikes 
are transient features. We will see that the initial Kasner point and the final Kasner point of 
a joint low/high velocity spike transition 7j are are related by the spike map (37'), which thus 
characterizes both Th\ and 7j transitions; we will argue that this property is not coincidental. 

In order to unveil the structures leading to joint low/high velocity spike transitions 7j we use 
information about timing. First, consider Fig. 4(d) and compare the trajectory of the spike surface 
z = with the trajectories associated with z 7^ of a neighborhood Uo of z = of <S|_o- The (dashed) 
lines of constant time show that the respective final Kasner points, A and A\, are approached 
'simultaneously' (which is in contrast to the approach to the initial Kasner state); specifically, 
let X(t, z) denote the family of curves representing the S\_ a transition; then \X(t, 0) — Aq\ and 
\X(t, 1) — Ai\ are of the same order as t — >• 00. Furthermore, for t » 1, the set {X(t, z) \ z € U } 
is approximated by the straight line Sa connecting the two Kasner points A and A\, which 
corresponds to the T/Vi transition that originates from A , see Fig. 5. 

Second, consider Fig. 4(a) or 4(b); let Y(t, z) denote the family of curves representing the %\\ 
transition. The convergence as t — > —00 to the initial Kasner point Bq is non-uniform; the 
(dashed) lines of constant time show that for t <C —1, the set {Y(t, z)\ z e R} is approximated by 
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the straight line Sb representing the T/Vi transition connecting the initial Kasner point B with 
the Kasner point B\. 

In the general Gi case the variable R\ does not vanish identically. Consider initial data generated 
by a small perturbation of <S L o initial data such that |i?i| 7^ is sufficiently small. (In accordance 
with the considerations of section 4.3 we assume symmetric initial data sets; in particular, R\ is 
an even function of z.) The associated solution S(t, z) shadows the S\_ Q transition closely; hence, 
for a large tA, {S(tA, z) \ z} rts Sa- However, S(t, z) cannot converge to the Kasner circle since R\ 
increases in a neighborhood of Sa, cf. Fig. 10. Eventually, despite the initial smallness of R\, the 
solution S(t,z) will be transported away from Sa by a 'frame transition'. 

A trajectory {S(t, z)\ t} with sufficiently large \z\ (i.e., sufficiently far from the spike surface) is ap- 
proximated by a BKL chain: While S(t, z) approaches the Kasner point Ai, R\ grows steadily and 
a Tr ± frame transition follows, i.e., S(t,z) ss Tr 1 for some time interval, which brings S(t,z) into 
the vicinity of the Kasner point B\, where R\ decreases rapidly; from there a Tr 3 frame transition 
ensues taking S(t, z) to the Kasner point Co where the trigger N\ induces a T/Vi curvature transi- 
tion that eventually takes S(t, z) to a Kasner point C\. For the trajectory {S(t, 0)\,t} representing 
the spike surface, the scenario is straightforward as well: S(t, 0) approaches the Kasner point A 
while Ri grows steadily; a Tr x frame transition follows, which brings S(t, 0) into the vicinity of 
the Kasner point Bo (where i?i decreases rapidly); a spike surface trajectory 7s follows and the 
solution approaches the Kasner point C\; in particular we observe that S(t, 0) rejoins S(t, z) at C\\ 
we refer to Fig. 12. The described scenario is supported by heuristic arguments invoking Fig. 10(b) 
and the symmetry of the initial data set, and by numerical simulations, exemplified in Figs. 12(a) 
and 12(b). 

Numerical experiments provide further information: For t > tA, the solution S(t, z) is transported 
from Sa to Sb, i.e., there is t B > tA such that {S(ts,z)\ z} rts S B , see Figs. 12(a) and 12(b). Note 
that i?i is decreasing rapidly as the solution approaches S B - The 'transition' from Sa to Sb occurs 
in a special manner: Curves of constant time, i.e., {S(t, z)\ z}, tA < t < t B , see the dashed lines 
of Fig. 12(b), are straight lines corresponding to segments of the paths of curvature transitions. 
The 'frame transition' that maps S(t,A,z) to S(t Bl z) is therefore the (nonlinear) superposition of 
two motions: The motion represented by the 'ladder' of constant time slices of Fig. 12(b) and the 
motion of spatial points along t = const lines. For t close to tA, (34a') yields 

Pl ~ f(u A ) + f(-u A )k\e^{t-t A ) z 2 ' 

where k A is a constant and u A is associated with the Kasner point A ; (51a) describes the motion 
of spatial points z along Sa for t close to tA- (Note the intimate connection of (51a) with a Bill 
spiky feature, see section 3.5.) For t close to t B , on the other hand, we have 

-w B (l-fc|e 4 ^(*^)z 2 ) _ -u A {l - kle^V-tzh 2 ) 
Pl ~ f(u B )+f(-u B )k 2 B e^t-t B ) z 2 - f(u A ) + f(- UA )k 2 B e^B(t-t B ) z 2 > ( 51b ) 

where k B is another constant and ub is associated with Bq, i.e., ub = u~ A (because ub arises from 
u A through a Tr 1 transition). Equation (51b) describes the motion of spatial points along Sb 
when t is close to ts- Interpolation between the two motions (51a) and (51b) (qualitatively) leads 
to the trajectories observed numerically, an example being the curve connecting Aq with B\ in 
Fig. 12(a). The remaining (and crucial) question concerns the evolution for t > tg. The numerical 
simulations show that the solution S(t, z) continues as an approximate high velocity transition 
T\\\ for t > ts (while R\ continues to decrease). In particular, we find that the entire family of 
trajectories S(t, z) rejoins at the final point C\ of the Th\ transition. 

Figs. 12(a) and 12(b) summarize the considerations and give depictions of a joint low/high velocity 
spike transition 7j - We reemphasize the property that the t = const curves are straight lines that 
coincide with the paths of Tn x curvature transitions. This 'timing' property shows that the spiky 
structure that forms under an (approximate) S\_ a is transported by a R\ frame rotation to the spiky 
structure appearing in (the middle of) an (approximate) %\\ and eventually 'un-forms' according 
to the 7hi evolution. 



(51a) 
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( a ) 7jo (b) Skeleton of 7j with t — const curves 



Figure 12: A joint low/high velocity spike transition 7] consists of a low velocity spike solution 
S\_ followed by a family of Tr 1 frame transitions and the 'second part' of a high velocity spike 
transition Th\. In (a) we depict the numerical evolution of four timelines associated with different 
values of z from z — (spike surface trajectory) to \z\ ^> 1 (BKL chain). Subfigure (b) shows the 
curves of constant time (dashed lines) for a 7j , which correspond to segments of the straight lines 
representing T/Vi transitions. 



A joint low/high velocity spike transition 7j induces a map between the initial and the final Kasner 
state of the transition, 

__ Th;, 
Slo I 1 — U- 1 7ri I i-u_ 



u_— 1 1 — u_ 

which, in terms of the standard Kasner parameter coincides with the spike map (37') of 7h\ tran- 
sitions. The connection between Th\ and 7j transitions is even stronger: Consider Fig. 14. On 
the one hand, the Kasner point O in sector (132) with 1/2 < u < 1 is the initial point of a joint 
low/high velocity spike transition 7j ; on the other hand, there is a Tr 1 frame transition emerging 
from O, which can be continued with a high velocity transition and another Tr 1 frame transition; 
the finite chain Tr 1 -Th\-Tr 1 yields 

u_ i— -> — i — > i— > = u + , (52b) 

m_ «_ 1 — u_ 

hence the final state of the chain Tr 1 -Th ~Tr 1 coincides with the final state of the joint low/high 
velocity spike transition 7j - This corroborates that 7j and 7h; are structures that are on an equal 
footing (through additional R\ frame rotations). 

A spike chain is an in general infinite concatenation of spike transitions and frame transitions; the 
'chain links' are high velocity spike transitions Th\ and joint low/high velocity spike transitions 
7j and Tr 1 and Tr 3 frame transitions. Examples of (short parts of) spike chains are given in 
Fig. 13. The role of spike chains for the dynamics of (generic) Gi models is suggested by numerical 
experiments: Asymptotically, solutions S(t,z) are approximate spike chains, i.e., solutions shadow 
spike chains (with an increasing degree of accuracy) . The role of spike chains for generic G 2 models 
is thus exactly analogous to the role of BKL chains for spatially homogeneous models. 

These considerations strongly suggest that the 'permanent' asymptotic features ('permanent spikes') 
of T 3 Gowdy models, which have been emphasized in the literature, are irrelevant for general G 2 
models and generic singularities: The permanent features of T 3 Gowdy are merely transient fea- 
tures in the generic context. Spatial structures form and 'un-form' recurrently; however, during 
these 'oscillations', the spike widths shrink to zero size. Oscillations are oscillations between Kasner 
states; the asymptotic dynamics of a spike timeline is represented by the spike trajectory (2 = 0) 
of a spike chain; the asymptotic dynamics of a generic spatial point is represented by a BKL chain 
(which corresponds to the trajectory on the 'boundary' \z\ — 00 of a spike chain). 



6 CONCLUDING REMARKS 



25 




Q 2 



(213) 



S 3 



y 




1123) 



(123) 



Qi 



Si 



Til 
(32l) 



Qi 



(132) 



(132) 



(a) Spike chain with two T\-\\ transitionsy 



(b) Spike chain with 7hi and Tj 



Figure 13: Examples of (parts of) spike chains. In (a) the continuous lines represent the spike 
timeline [z = trajectory) of 7h~Tr 3 -7r 1 - T\\\—Tr 3 -Tr x — Tr 3 ; the dashed lines are Tr 1 -Tr 3 -Tr 1 ■ 
In (b) we have Th;-Tr 3 -Tj (where 7j consists of three pieces: Slo-Tr^Thi) ■ 

We conclude with two remarks. First, consider Fig. 14 and recall that the final state of a low/high 
velocity spike transition 7j and the chain Tr^ThtTr^ coincide. It is highly plausible that 'mixed' 
(or 'double') transitions exist, in close analogy with, e.g., double frame transitions (involving R\ and 
i?3 at the same time) and mixed curvature/frame transitions (involving iVi and frame variables at 
the same time), see [48]. This further stress the close relationship between 7j and Th\ and identifies 
7j and Tr^Th-Tr^ as the limiting cases (boundary) of an entire family of transitions. However, 
by analogy with the case of double frame and mixed curvature/frame transitions we expect that, 
asymptotically and generically, merely the single spike transitions, 7j and Th\, will occur in a spike 
chain. Second, we note that in the spatially homogeneous cases Vl_i/ 9 , VIII, and IX, solutions 
cannot converge to Kasner states (other than the Taub points — in the LRS case); all solutions 
are oscillatory. We expect (almost) the same to be true in the general G2 case: While in the OT 
case, solutions converge to generalized Kasner states, we expect that (almost) all solutions are 
oscillatory in the general G 2 case. 

6 Concluding remarks 

In this paper we have discussed the asymptotic dynamics of G 2 models toward a singularity. 
We have given heuristic arguments and strong numerical evidence showing that the well-known 
permanent spatial structures (spikes) arising in special models such as the T 3 Gowdy models 
are absent in general G2 models: General G2 models possess oscillatory singularities, which are 
characterized by the recurring formation and 'un-formation' of spiky features whose width rapidly 
shrinks towards the singularity; the asymptotic dynamics is represented by infinite spike chains 
and thus associated with sequences of Kasner epochs. 

The world line of the spike surface is characterized by an oscillation between Kasner epochs induced 
by the map that in terms of the gauge-invariant Kasner parameter u is given by 




u-2 



(u-l)- 1 -! 




u e [3, 00) , 
uG [2,3], 
u €[3/2,2], 
ue [1,3/2], 



(53) 
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Figure 14: A Kasner point O in sector (132) with 1/2 < u < 1 is the initial point of a joint low/high 
velocity spike transition 7j (where the continuous line represents the timeline of the spike surface) ; 
however, there is also a Tr 1 frame transition emerging from O, which can be continued with a Tn x 
curvature transition and another Tr x frame transition. The two trajectories meet at the same final 
Kasner point. 

which is in contrast to the usual Kasner map 



In a subsequent paper [69] we will we analyze statistical properties of the Kasner oscillations of 
the spike world line by investigating the map (53) in detail. 

Although we expect that non-local (i.e., non-BKL) dynamics will be asymptotically confined to 
spike surfaces, which is due to the rapidly shrinking spike widths, and that these spike surfaces 
constitute a set of measure zero of all spatial points, an understanding of recurring spikes is still 
crucial for our understanding of generic spacelike singularities: 

• The key solutions that describe BKL and non-BKL behavior are intimately related with 
each other through solution generating algorithms, which create a hierarchical order for 
these solution; this strongly hints at the existence of hidden symmetries, which are thus 
a common theme in our quest to understand generic spacelike singularities. 18 The case for 
hidden symmetries is further strengthened by the observation that the 'spatial boundary' (i.e., 
\z\ 00) of the spike chains are the BKL chains. Since the spike surface trajectories (non- 
BKL dynamics) meet the trajectories of the other spatial points (BKL dynamics) periodically, 
there is a close relation between the map (53) and the Kasner map (54). Accordingly, the BKL 
scenario seems to be part of a greater picture which is needed for a complete understanding 
of generic singularities. 

• Proofs about generic spacelike singularities must take into account recurring spikes, since 
estimates (e.g., of spatial derivatives) are heavily affected. 

18 The solution generating techniques used to obtain the spike solutions in a hierarchical manner involve the OT 
line element in an Iwasawa frame, i.e. (47), and in the area time gauge. The solutions are obtained by alternately 
applying the so-called Gowdy-to-Ernst (GE) transformation and a special frame rotation (FR) (which is obtained 
by a rotation of the symmetry coordinates x and y by — 7r/2). Starting with the Kasner solution in diagonal form 
and performing a GE transformation yields the Kasner solution in non-trivial OT form; a FR yields the frame 
rotated Kasner solution associated with Tr 3 ; subsequent application of a GE transformation gives the Bianchi type 
II vacuum solution in a Fermi frame; applying the FR again yields a 'spikily rotating' Bianchi type II vacuum 
solution, known as a false spike solution; acting on this with the GE transformation results in the explicit spike 
solution; for details, see [53]. 
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• Spikes are associated with the zeros of certain functions (e.g., the variable Ni). Special initial 
data (like the symmetric data considered in this paper) can fix the location of a spike; however, 
in general, spikes need not be present initially but will form when a function goes through zero 
and move in accordance with the evolution of the functions and its zero(s). Unfortunately, 
this does not answer the question what the physical reasons for spike formation are. Our 
lack of knowledge and intuition for spike formation prevents us from even an educated guess 
about how many spikes form generically. At the moment we cannot exclude that a dense set 
of spikes can form toward generic spacelike singularities. BKL asymptotic dynamics would 
still be generic, but the BKL scenario would certainly obtain an unexpectedly significant 
non-BKL counterpart. 

Let us turn to numerical issues. At present, numerical accuracy regarding the simulation of spikes 
can be categorized into three levels: 

• At the crudest level of numerical accuracy there are not enough grid points to simulate spikes 
correctly; moreover, spikes are artificially produced and annihilated. 

• At the intermediate level of numerical accuracy there are enough grid points to accurately 
describe one spike transition correctly, but numerical convergence is not achieved for simula- 
tions that are supposed to follow a concatenation of spike transitions. The first indication of 
non-convergence is the difference in the timing of transitions. The intermediate level yields 
quantitatively correct results for simulations of one spike transition, and qualitatively correct 
results for simulations of short spike chains. 

• Finally, at the highest level of numerical accuracy there are enough grid points to obtain 
numerical convergence for longer simulations that cover short spike chains. But even these 
simulations are currently limited to two or three spike transitions. We are confident that 
improvement can come from the heuristic insights gained in this paper. 

With the benefit of hindsight we are in a position to reproduce and assess what was actually 
achieved in the numerical simulations of different papers. Early numerical work on 'spiky features' 
has to be considered as having crude numerical accuracy, thus belonging to the first category, 
but it nevertheless managed to tie spiky features to analytic conditions in special models such as 
the T 3 Gowdy models, as exemplified in [25]. Hern [31] and Berger et al (see [10] and references 
therein) also studied general Gi models with T 3 topology numerically, but again, even though Hern 
used mesh refinement, numerical accuracy was not sufficient to resolve several spike transitions 
correctly. 19 It was not until 2003 that the first qualitatively correct spike simulation was done by 
Garfinkle and Weaver [36] in the case of T 3 Gowdy models, described by means of P T and the 
speed v (see Appendix A. 4 and Fig. 8 for a description in terms of projections on (Xi, £2, 
space). The first qualitatively correct spike simulation in the general Gi case was first done in 
2005 by Andcrsson et al [50]. This qualitative picture was confirmed quantitatively and linked to 
the Lim's explicit solutions [53] in 2009 by Lim et al [54]. The key idea to simulate several spike 
transitions accurately was to use the explicit spike solutions to design a grid that zoomed in on a 
single recurring spike. Despite this, only a few high (but no low) velocity spike transitions were 
followed. 

The work [54] makes it clear what a formidable numerical challenge recurring spikes pose. It is 
important to note that all quantitatively correct spike simulations so far involve recurring spikes 
that are forced by special initial conditions to be fixed in space (non-moving spikes). The spike 
surface thus obtains a fixed value of the spatial coordinate, e.g., zspiko = 0. In this paper we have 
pushed the current numerical frontier in the context of the recurrence of ' 'non-moving spikes' by 
numerically correct simulations of a chain of spike transitions that also involve low velocity spike 
solutions. The main challenge in this context is that low velocity spike solutions require numerical 
accuracy at super horizon scales. 

In general, recurring spikes will not be spatially fixed, instead they will move, i.e., zspikc = zspike(i)- 
These 1 moving spikes' are not directly described by the explicit spike solutions (30). However, it 

19 Hcrn also simulated spikes in models with a single spacclike Killing vector [31]. 
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is possible that moving spikes 'asymptotically freeze' in space which implies that the description 
of the recurrence of these spikes in terms of spike chains straightforwardly applies. 20 The question 
of asymptotic freezing is closely related to another unexplored area: Present numerical techniques 
are not sufficient to accurately describe the creation and the possible annihilation of spikes. This 
makes it difficult to guess whether the number of spikes that form remains finite, or whether, in 
general, an infinite number of spikes, possibly a dense set, is created in the asymptotic regime. It 
is doubtful whether a sound numerical investigation of this question is possible unless the numerics 
is supported by analytical insights. 

There are many examples in general relativity that show that one must go beyond a given context 
in order to understand it, especially as regards asymptotics. This is particularly relevant as regards 
topology. For instance, Bianchi type I and II models are essential building blocks for the under- 
standing of the singularities of the more general Bianchi models; it is important to note that this 
is irrespective of the fact that these Bianchi models have spatial topologies that arc completely 
different from those of Bianchi types I and II. Yet another example is provided by the present 
paper: Understanding the singularities of G2 models requires to go beyond T 3 topology, since the 
essential building blocks, the explicit high and low velocity spike solutions do not admit such a 
topology. Special models like these provide local building blocks (where local may refer to the 
particle horizon scale) for the asymptotic description of more general models, or even generic ones. 

The relationship between local and global issues deserves special attention. It is important to realize 
that the primary importance of the G2 models as regards generic singularities is not the models 
themselves, but the crucial property that the G2 equations are those of an invariant boundary 
subset of the general conformally Hubble-normalized state space (associated with the Einstein 
equations without symmetries), which we call the partially local G 2 boundary. 21 This subset is 
obtained by setting certain frame variables to zero to achieve d±(*) = 9 2 (*) = in the field 
equations, where d\ and 82 are Hubble-normalized spatial frame derivatives; by this, the equation 
induced on the partially local G2 boundary are the same as those for the G2 models, where, 
however, the constants of integration may now depend on the spatial coordinates x and y. The 
situations is thus completely analogous to the relationship between the local boundary and the 
spatially homogeneous models. Since what survives of the G2 models in the context of the general 
conformally Hubble-normalized state space is merely the equations (on a boundary subset), any 
global considerations, which may be of interest for the G2 models themselves, are probably not 
particularly relevant in the general context (for general models there typically are other global 
considerations that are relevant). As regards singularities it is the local results, concerning domains 
of dependence in the vicinity of the singularity, that are of importance for, e.g., cosmic censorship 
in the generic case for models without symmetries. For further discussion on topological issues in 
the case of G2 models, see Appendix A.l. 

Going beyond the G2 assumption and looking at models with fewer isometries does not only shift 
the attention to the partially local G2 boundary, it leads to completely new challenges as well. Spike 
surfaces, within the context of the partially local G2 boundary, are no longer planes, and they can 
intersect in curves that lead to different spike dynamics. Similarly, if one has no symmetries 
at all, spike curves may intersect at points, which again may lead to new spike dynamics. At 
present it is not known whether such intersections persist or recur, or if they are transient and 
hence irrelevant asymptotically. 22 The above issues illustrate again that spike dynamics poses a 
formidable numerical and analytical challenge. It is clear that we are only beginning to understand 
the rich structure of generic singularities, and the underlying physical reasons for their existence 
and characteristics. 

20 There exists heuristic evidence for asymptotically frozen spikes. In the billiard picture, a 'dominant wall' that 
has a hole in it at an isolated value of a spatial coordinate corresponds to a spatially frozen spike [49] . 

21 The partially local G2 boundary, previously referred to as a the partially silent boundary subset, yields, e.g., 
spiky behavior that is still associated with asymptotic silence, but also solutions with singularities that break 
asymptotic silence, which are characterized by the existence of directions in which particle horizons extend to 
infinity, see [64]. 

22 Weak numerical evidence suggests that intersections only occur momentarily, and hence that the BKL picture 
in combination with G2 spike oscillations might capture all essential features of generic spacelike singularities. 
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A Appendices 

A.l G 2 models and topology 

A G 2 model is a spacetime admitting a two-dimensional Abelian isometry group. In the literature 
the spacetime is often assumed to be globally hyperbolic and the action to be effective on spacetime 
Cauchy hypersurfaces with two-dimensional principal orbits. A substantial part of the rigorous 
results concerns G 2 models with a compact Lie group which is identified with T 2 = U(l) x U(l); 
for this reason these models are occasionally referred to as T 2 -symmetric spacetimes'. By assuming 
compactness of the Cauchy hypersurfaces it follows that the topology has to be T 3 , S 2 x S 1 , S 3 , 
or that of one of the Lens spaces L(p; q) (which are quotient spaces of S 3 ), see [70]. In the cases 
of S 2 x S 1 , S 3 , and L(jp; q), the group action must have fixed points (which correspond to zeros of 
the Killing vector fields), and hence the two 'twist constants' 

c A = vol abcd X?X*V c Xi (A = 1,2) (55) 

of the two commuting Killing vector fields Xi and X 2 vanish. (The twist constants are spacetime 
constants, i.e., dcA — 0, if the Einstein vacuum equations are imposed [68, 70].) In other words, in 
these cases the action is necessarily 'orthogonally transitive' (OT) [7], i.e., the 2-spaces orthogonal 
to the group orbits are surface forming. 

The spatially compact topology that is of particular interest in the context of generic singularities 
is T 3 = S 1 x T 2 (and its covering R x T 2 ), because in this case the twist constants do not vanish 
in general, i.e., in the T 3 case, c\ = c 2 = is a restriction. The (standard) T 3 case arises from 
the R x T 2 case by identifying {0} x T 2 and {2n} x T 2 by means of periodic boundary conditions 
for the metric and the extrinsic curvature. (The corresponding OT models are the T 3 Gowdy 
models [24].) It is also possible to consider non-trivial torus bundles over S , which correspond 
to identifying {0} x T 2 and {2tt} x T 2 through non-trivial SL(2,Z) transformations. In this case, 
the action of T 2 is global on the covering space K x T 2 but merely local on T 3 (i.e., there do not 
exist global Killing vector fields), see, e.g., [72], 23 and the metric and the extrinsic curvature satisfy 
non-trivial boundary conditions. 

It is possible to w.l.o.g. assume the vanishing of one of the twist constants; namely, if c\ 7^ and 
c 2 7^ 0, there are linear combinations X[ and X' 2 of X\ and X 2 , which are again Killing vector fields, 
such that the associated twist constants are c x = and c y 7^ 0. We introduce local coordinates 24 
such that X[ = d x and X 2 = d y . Then the metric takes the form (1'), which coincides with the 
form in [30, 57] by setting U = -b 1 , logR = -b 1 - b 2 ,v = -b 1 - 6 3 , a = N 2 exp(2& 3 ); the shift is 
set to zero, see also [58, 59]; finally, A = — n\, G\ = —n 2 , G 2 = —T13. In the spatially compact case, 
these functions are periodic in z (or, in the case of a non-trivial torus bundle, satisfy boundary 
conditions derived from the particular SL(2, Z) identification). In terms of the metric functions, 
cf. (1), and the area density R, the twist constants are represented by 

c x = , (56a) 
c y = R 3 e 2bl+b3 N- 1 d x on 3 = e^-^W N -i dx0Tl3 . ( 56b ) 

The vanishing of the first twist constant is equivalent to the identity 

-d x on 2 + nid x ons = d x on 2 - n^d Xa ni = , (56c) 

23 Thesc considerations arc also relevant in connection with the specialization of G2 models to spatially homoge- 
neous Bianchi models [72]. 

24 Unless X\, X2 and X[, X' 2 are related by an 5L(2,Z) transformation, x and y are not standard coordinates on 
T 2 ; however, this does not affect the equations; see [70]. 
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i.e., i?2 = 0, see eq. (58f) below. In area time gauge (3) the metric takes the form (4), where a 
satisfies the equation 

d t a = e 2 » 2 - 2b3 a 2 c 2 y , (57) 
see [57], and c y becomes c y = — e~ 2f> +2b oT x l 2 dtn 3 . 

The OT case is characterized by vanishing twist constants; through (56) this leads to functions n 2 
and n 3 that are independent of time. Using the coordinate freedom x i-> x + fi(z), y i-> y + f 2 (z) 
we are able to set n 2 = n 3 = without loss of generality. Likewise, we obtain that a is independent 
of time from (57). Using the remaining coordinate freedom z ^ fz{z) and redefining b 3 we may 
thus set a = 1, which corresponds to a lapse function (5). Note that setting a = 1 is impossible 
in the general case (with c y ^ 0) as long as the area time gauge is enforced; on the other hand, 
abandoning (3) makes the 'conformal gauge' (or: 'null cone gauge' [61]) a = 1 possible [57]. 



A. 2 Metric variables and Hubble-normalized variables in the G2 case 

In the G 2 case, the general relations of Appendix A of [48] reduce to the following relationship 
between the metric variables b 1 , b 2 , 6 3 , n 1; n 2 , n 3 , and the Hubble variable H and the conformally 



Hubble-normalized state space variables: 

E 3 = H- 1 exp(b 3 ) , H = -lN- 1 d x o(b 1 +b 2 + 6 3 ) , (58a) 

r = -H- 2 exp(b 3 )d z H, U = -r + {NH)- 1 cxp(b 3 )d z N , (58b) 

A = r+ iff- 1 cxp(6 3 )^(fe 1 +6 2 ), q = -1- N^H^d^H. (58c) 

For simplicity we have set E3 = E^, A = As, r = r^, and U = U 3 . These quantities refer to an 
Iwasawa frame; functions depend on x° and z alone. Furthermore, 

Y, a = -1- {NH)- 1 d x o{b a ) {a = 1,2, 3), (58d) 

i?i = -£ 23 = -i(Aff)- 1 cxp(fe 3 - b 2 )d x o(n 3 ) , (58e) 

i? 2 = E 31 -i(Ai/)- 1 cxp(6 3 -6 1 )[a ;c o(n 2 )-n3a x o(n 1 )] , (58f) 

R 3 = -E 12 = -liNH)- 1 cxp(fe 2 - b^d^im) , (58g) 

Ax =i/- 1 exp(& 2 +6 3 -6 1 )9 z (n 1 ), N 12 = cxp^d^b 1 - b 2 ) . (58h) 



Adapting the spatial coordinates to achieve the vanishing of the first twist constant leads to 
d x on2 — n 3 d x ani = 0, see (56c), and thus i? 2 = 0. In area time gauge, the sum b 1 + b 2 does not 
depend on z, hence A = r in that case. The lapse N is unspecified in general; for OT models in 
area time gauge, however, we have N = — cxp(— b 1 — b 2 — b 3 ), see (5). 



A. 3 Vacuum field equations in the G 2 case 

The vacuum field equations for the conformally Hubble-normalized variables in the G 2 case are 
conveniently divided into decoupled equations and a coupled system of evolution equations and 
constraints. 

The decoupled equations (for the metric variables and the Hubble scalar) read 

d b a = -(1 + S Q ) , (a = 1,2, 3) d ni = -2i? 3 cxp(6 2 - b 1 ) , (59a) 

d n 2 = n 3 d ni , d n 3 = — 2i?x exp(6 3 — b 2 ) , (59b) 

d H = -(l + q)H, (59c) 

dsib 1 - b 2 ) = 2A 12 , dsib 1 +b 2 ) = 2(A-r), (59d) 

d 3 m = mexpib 1 -b 2 ), d 3 N = (U + r)N, (59e) 

d 3 H = -rH , (59f) 
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where the conformal frame vectors are 

do = {NH)- x d xQ , d 3 = E 3 d z . (60) 

The coupled evolution equations for the conformally Hubble-normalized variables are 

d Q E 3 = (g - E 3 )£ 3 , ( 61a ) 

d ^i = -(2 - g)Ej + 2i? 2 _ 3 5i _ i (a3 + A + 3Nl2 )(fi + 2r ) - \(U 2 - 2r 2 ) , (61b) 
d £ 2 = -(2 - g)E 2 - 2R\ + 2R\ - 5 S 2 -\(d z + A- 3N 12 ){U + 2r) - |(C7 2 - 2r 2 ) , (61c) 

O E 3 = -(2 - g)E 3 - 2i?f - 3 S 3 + §(d 3 + A)(U + 2r) + §(l> 2 - 2r 2 ) , (61d) 

a i?i = -(2-g + S 2 -S 3 )ii 1 , (61c) 

a i? 3 = -(2 - g + Si - E 2 )i? 3 + 3 S 12 - (17 + 2r) , (61f) 

a A = (g - E 3 )A + i(d 3 + J7)(2g + E 3 ), (61g) 

doA^i = {q + 2Si)JVi - 2(d 3 + L> + 2iV 12 )i? 3 , (61h) 

a iV 12 = (g - E 3 )/V 12 - i(0 3 + L>)(S! - E 2 ) . (61i) 

The coupled constraint equations take the form 

= 1 - E 2 - Q k - i(2d 3 -4A + r)r, (62a) 

= (d 3 -3A + 2r + N 12 )R 1 , (62b) 

= 2r + (d 3 - 3A + 2r)E 3 + (E x - E 2 )/V 12 + . (62c) 

In these equations we have used the abbreviations 

V k = ± (N 2 + AN 2 2 ) - | (2d 3 - 3A)A , (63a) 

3 51 = | (Ni + N 2 2 ) - 2iV 12 A - 18 3 (A - 3N 12 ) , (63b) 

3 5 2 = -±(iV 2 - 2iV 2 2 ) + 2iV 12 A - \8 3 {A + 3N 12 ) , (63c) 

3 5 3 = -i(iV 2 +4iV 2 2 ) + |a 3J 4, (63d) 
3 5 12 = -1(83- 2N 12 - 2A)N 1 , (63e) 

q = 2Y? - I [a 3 + l> - 2(A - r)] (£7 + r) , (63f) 

E 2 = i (E 2 + E 2 + E 3 + 277? + 2i? 2 ) . (63g) 



The conformally Hubble- normalized equations are not the most useful equations in the G 2 context. 
Instead one can adapt to the special structure of these models and use, e.g., conformally area- 
expansion-normalized variables, see [61, 51]. However, since we are interested in how the G 2 
case fits into the general context without symmetries we have focused on the conformal Hubble- 
normalized approach in this paper. 

A. 4 OT models: Asymptotic velocity dominance 

In this section we give a brief overview of existing rigorous results for T 3 Gowdy vacuum models. 
Since these results are intimately connected with the concept of (asymptotic) velocity, we begin 
by introducing this particular diagnostic tool. The metric (47) in area time gauge is sometimes 
written in the form 

ds 2 = - e (*- A )/ 2 (- e - 2 *^ 2 + dz 2 ) + e p - 1 (dx + Qdy) 2 + e-^dy 2 , (64) 

see, e.g., equation (15) in [37], where (t,x,y,z) corresponds to (r,a,8,ff). Obviously, (64) is of the 
form (47) with 2b 1 = -P + t, 2b 2 = P + t, 2b 3 = A/2 - t/2, m = Q. There exists a number of 
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diagnostic tools connected with (64), which were originally developed for T 3 Gowdy models but 
are useful for all OT models in area time gauge. In particular one can define kinetic and potential 
energy densities according to [67] 

/C = /C(M) = (P 4 ) 2 + (e P Q,t) 2 , V = T(t,z) = (e- t P, z ) 2 + (e p - t Q, z ) 2 , (65) 

and a velocity v as the square root of IC; more appropriately, the velocity should be referred to as a 
speed (since it is non- negative by definition). The quantities IC,V, and v are intimately connected 
with the Hubble-normalized variables. In [71] it was noted that the velocity v is 



v = v(t, z) = VlC = z TyT " = V3 V , 4 (66) 

i - 2^3 -L - 2^3 



At the same time we obtain 



_ 2(1 - P. t ) . 2(1 + P t ) _ -l+JC + T 

Pl = 3 + K + V> P2= 3 + IC+V> P3= 3 + /C + V ' (67a) 

^ = 6 3^TP> * = 12 jm^' ^ = " 6 3^' <™» 
from which we infer that 

# + 8 + % = 1 - 8 (3 + K + V)» E " 1 - 12 (3TWF' (67C) 

A considerable part of the analytic rigorous work on spacclike singularities has been influenced by 
the work by Eardley et al [73], which has resulted in a focus on asymptotic velocities. The term 
asymptotic velocity refers to the limit 

Voo — v oc (z) = lim v(t, z) = lim y/fC(t, z) = \1u + 1| . (68) 

t— >oo t— >oo 

In the context of the T 3 Gowdy models, the existence of this limit has been proved (but we 
emphasize that this limit does in general not exist for G 2 models without the OT assumption, see 
section 5). Specifically for each smooth solution (P, Q), is an upper semicontinuous function of 
z [67]. Furthermore, e 2P Q 2 t and V converge to zero as t — > oo and P 2 — > (t —¥ oo). (This holds 
not merely pointwise in z but uniformly on the shrinking particle horizons [z — e~ l , z+e 1 } associated 
with a z = const timeline; we refer to our discussion of particle horizons in section 4. Evidently, the 
spatial topology is irrelevant in this context.) Following Moncrief and coworkers, the asymptotic 
velocity Vca has been used extensively as a diagnostic tool see, e.g., [26, 34, 37, 67, 74, 75]. To 
resolve the ambiguity in the limit of P t , following [33, 34] one may define the asymptotic velocity 
component 

k = lim P t = lim - 1 . 2 - = 2u + 1 , (69) 

i.e., \k\ — Voo; it follows that < k < 1 corresponds to — \ < u < 0, i.e., to the stable sector (312). 

The convergence statements are commonly subsumed under the term 'asymptotic velocity domi- 
nance', and the existence of the limit (68) is referred to as convergence of solutions to an asymp- 
totically velocity dominated state. As a consequence of (67), however, the convergence statements 
simply mean that each solution converges to a generalized Kasner metric with a semi-continuous 
Kasner parameter u — u(z). 

The convergence of solutions to the Kasner circle K° is understood in quite some detail in the T 3 
Gowdy case. Building on [26] and the results of Rendall and coworkers [34, 35], Ringstrom proved 
that solutions (P, Q) with < k < 1 (which implies k = v^) possess asymptotic expansions of the 
form 

P(t, z) = Voo {z)t + 4>{z) + u(t, z) , Q(t, z) = q(z) + e- 2v ~ {z)t {ij{z) + w{t, z)) , (70) 
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where u and w and its derivatives converge to zero exponentially, see Prop. 1.5 in [67] for the 
precise statement. 25 Note that the expansion (70) represents the convergence of solutions to the 
stable sector of the Kasner circle. In [75, 76, 77] criteria on initial data have been derived that 
guarantee the assumption < k < 1 on the asymptotic velocity. By these theorems, the heuristic 
reasoning of section 4.2 that solutions that are sufficiently close to a generalized Kasner metric on 
the stable sector converge to the stable sector is made rigorous in the context of T 3 Gowdy models. 

Besides (70) there exist solutions with asymptotic expansions of the type (70) where, however, 
the asymptotic velocity is discontinuous. These are the solutions that exhibit a 'true' spike in 
the nomenclature of [37]. 26 Alternatively, these solutions are said to possess a 'permanent' spike. 
The explicit low velocity spike solutions (38) fall into this category, while the explicit high velocity 
spike solution (37) possess a constant (and thus continuous) asymptotic velocity. However, we 
strongly emphasize that the 'permanence' of low velocity spike solutions, i.e., convergence (to a 
discontinuous limit), is true in the OT context only. As seen in section 5, spikes are not permanent 
but transient and recurring features in the general Gi case. 

We note that the usefulness of the (asymptotic) velocity in the characterization of the behavior of 
OT models is indisputable since it captures essential aspects of the dynamics of solutions; however, 
it should not come as a surprise that v fails to capture all aspects of the dynamics. Furthermore, 
the role of the velocity is diminished when we go beyond the OT case where the asymptotic 
velocity does not even exist since v does not converge, see section 5. However, we obtain a 
clearer picture when we do not restrict ourselves to one particular degree of freedom but instead 
consider projections of the Hubble-normalized state vector onto the 2-dimensional space spanned 
by (Si, £2, £3)) as is done, e.g., in the present paper. 

Finally, the issue of 'high' and 'low' velocities deserves a comment. In the Gowdy literature 
the question of whether a velocity is high or low usually refers to the asymptotic velocity Woo. 
However, if the limit (68) does not exist, as for the general G 2 models, this classification becomes 
meaningless. In addition, the explicit spike solutions (30) do not easily fall into the 'high/low' 
asymptotic velocity scheme: For the explicit high velocity solutions the limit is an arbitrary 
positive number; for the explicit low velocity solutions the limit is in the interval (1,2) for the 
spike timelines and in the interval (0, 1) for the z ^ timelines. In the context of the explicit spike 
solutions it is thus preferable to classify solutions in terms of their limit t — > —00; the corresponding 
(anti) asymptotic velocity (as t — > —00) is in the interval (2, 3) for the explicit low velocity spike 
solutions and in the interval (3, 00) for the explicit high velocity spike solutions, thus motivating the 
present nomenclature for the spike solutions. With hindsight, and with issues concerning generic 
singularities and cosmic censorship as ultimate goals, we find that there is room for improvement 
concerning the terminology that is used in much of the literature. 



A. 5 Evolution of spike widths 

In this section of the appendix we show the result given in equation (50). Consider equation (61f). 
In area time gauge, since (58c) implies that A = r, we have 

Af^dtRs = -(2 - q + E x - S 2 )i? 3 - \ AsM + JV1JV12 - \N\V , (71) 

where U = d$ log (—AC) and where the Hubble-normalized lapse N is given by 

U = NH = ~2^%> (?2) 

which follows from b 1 + b 2 = t and (58d) by using that S3 = — Si — £2- Making use of (63f) and 
the constraint (62a) equation (61d) turns into 

AA- 1 a t S 3 = (2- g )(2-S 3 )-2i? 2 . (73) 

25 The assumption < fc(zo) < 1 for some zq is sufficient to obtain smoothness of Voo and uniform expansions of 
the type (70) in a neighborhood of zq. 

26 In this context, the term 'true' refers to the non-uniform convergence of scalar functions, e.g., the Hubble- 
normalized Kretschmann scalar, to a discontinuous limit. 'False' spikes (which we refrain from discussing) possess 
a discontinuous k, but v^o and curvature scalars arc continuous. 
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Therefore, in the OT case, where Ri = 0, we hnd that 

N^dtN = (2 - q)N . (74) 

Accordingly, equation (71) yields 

M- 2 d t (MRs) = -(Si - S 2 )i? 3 - \ E 3 d z N t + iViJVia - 9 Z log(-JV) . (75) 

Furthermore, from (61a) and (74) we obtain 

d t (AfE 3 ) = -NE 3 , (76) 
which entails that AfE 3 = — fce~* for some k — k(z); we set n = fc(0). 

Consider an OT solution S that is an approximate OT spike chain with spike surface z — and 
symmetric functions, i.e., Si, S 2 , S 3 , and R 3 are even and Ni, A 12 are odd. We evaluate (75) at 
the spike surface z = to obtain 

d z N x \ z=Q = -2(AA^ 3 )- 1 | z=0 (Ar~ 1 d t (AfR 3 ) + (Si - S 2 )AO? 3 ) \ z=0 ■ (77) 

Insertion into (49) yields an expression for the spike coordinate width c z according to 

c z = (A/'-E 3 )~ 1 | Q x function (spike surface orbit z = 0) . (78) 

Note that the spike surface orbits at z — are independent of the width of the (approximate) 
spike transition. For an explicit spike solution with width c and time offset —t, i.e., (30) where z 
is replaced by cz and t by t — i, we have ME 3 = —cr 1 e~ lyt ~ t \ which follows from (58a) and (42), 
while the function of the spike surface orbits yields — e~(*~*). 

Suppose that the approximate OT spike chain is approximated in some time interval by an exact 
spike transition with parameter uq, (inverse) width cq, and time offset — i ; hence, by (78), 



Co 



= -K-Ve-C-'o) . (79) 



At a (much) later time, when the solution is approximated by the subsequent %k of the OT chain; 
the Kasner parameter is u\ = uq — 2, the (inverse) width is c\ and the time offset is —t\ (where 
ti 3> to). According to (78) and (79) we have 

ci = -k-V e-C-' 1 ' = coe*-* 1 ^-**-* 1 ) = coe 1 ^ , (80) 

i.e., the (inverse) coordinate width increases by a factor of e* 1_t °. This establishes (50). 
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